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We  have  simulated  the  wind  measurements  which  might  be  obtained  from  a  low  power  doppler  wind 
lidar  (DWL)  sensor.  Wind  profiles  obtained  from  such  an  instrument  will  not  be  as  complete  as  those 
expected  from  planned  higher  energy  sensors  (e.g.  LAWS),  with  measurements  only  from  citrus 
cloud,  the  tops  of  other  clouds,  and  from  marine  boundary  layer  aerosol.  The  potential  impact  of  such 
a  scaled  down  DWL  system  has  been  as.sessed  by  inserting  the  simulated  measurements  into  the  Air 
Force  Geophysics  Laboratory  (GL)*  Global  Data  Assimilation  System  (GDAS).  Results  from  the 
observing  system  simulation  experiment  (OSSE)  indicate  that  significant  improvement  in  the  analysis 
and  forecast  of  geopotential  height,  horizontal  wind,  and  relative  humidity  may  be  expected  in  the 
.southern  hemisphere  extratropics.  Compari.son  of  impacts  with  those  obtained  in  a  previously 
conducted  OSSE  which  assumed  a  full  blown  DWL  instrument  (WINDSAT)  indicate  that  the  low 
energy  lidar  yields  improvements  which  are  1/2  to  3/4  the  magnitude  of  improvements  with 
WINDSAT.  Thus,  if  the  development  cost  of  a  .simple  lidar  is  found  to  be  substantially  less  than  a 
LAWS-type  sen.sor,  such  an  instrument  may  warrant  consideration  in  the  planning  of  future  space- 
based  observing  systems. 
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1.  INTRODUCTION 

Tliis  report  summarizes  our  efforts  to  assess  the  usefulness  of  a  doppler  wind 
lidar  (DWL)  instrument  for  numerical  weather  prediction.  Recently,  a  great  deal 
of  attention  has  been  focu.sed  on  the  remote  sensing  of  winds  from  space  (Curran 
and  others,  1988).  The  propo.sed  Laser  Atmospheric  Wind  Sounder  (LAWS) 
instrument  would  allow  the  retrieval  of  wind  profiles  down  to  cloud  tops  globally. 
However,  LAWS  has  very  significant  energy  requirements.  In  this  study  we 
consider  the  impact  of  a  low  energy  DWL  instrument.  Such  an  instrument  will 
not  be  able  to  retrieve  winds  in  clear  air  but  will  obtain  winds  from  cirrus  layers, 
the  marine  planetary  boundary  layer  and  from  cloud  tops.  We  evaluated  the 
impact  of  such  an  instrument  by  simulating  its  data  coverage  and  accuracy,  and 
performing  an  ob.serving  system  simulation  experiment  (OSSE)  in  which  all  data 
are  simulated  and  inserted  into  the  Air  Force  Geophysics  Laboratory  (GL)  Global 
Data  Assimilation  System  (GDAS).  Resulting  analyses  and  forecasts  are  compared 
to  those  in  a  control  scenario  in  which  lidar  data  is  withheld  to  determine  impact 
of  the  DWL  data. 

This  report  summarizes  the  results  of  this  experiment.  Following  the  introduc¬ 
tion,  a  short  background  of  the  problem  is  discussed  in  section  2.  This  includes  a 
brief  historical  perspective,  recent  developments,  as  well  as  the  motivation  for  the 
present  study.  The  general  technical  approach  is  discussed  in  section  3,  followed 
by  the  specific  design  of  the  OSSE  in  section  4.  The  simulation  methodology  for 
the  DWL  data  is  covered  in  detail  in  section  5.  Section  6  contains  a  discussion 
of  the  data  impacts  In  the  analy.scs  and  forecasts.  Finally,  section  7  contains 
summary  conclusions  and  sugge.stions  for  future  work. 

2.  BACKGROUND 

Since  the  advent  of  meteorological  satellites  a  quarter  of  a  century  ago,  our 
space-based  perspective  of  the  earth’s  atmosphere-occan  system  has  been  acquired 
almost  exclusively  employing  passive  imagers  and  sounders.  From  the  earliest  sin¬ 
gle  channel  cloud  imagers  to  contemporary  multispectral  .scanners  and  sounders, 
passive  sen.sor  systems  have  undergone  continuous  evolution  with  each  succeeding 
generation  of  instrument  resulting  in  higher  spectral  resolution,  narrower  fields- 
of-view,  and  improved  detector  .sensitivity.  These  advances  have  been  manifested 
in  more  accurate  and  spatially  resolved  measurements  of  desired  meteorological 
fields  such  as  those  of  cloud,  temperature,  moisture  and  winds.  To  date,  passive 
systems  have  been  the  key  components  of  the  deployed  meteorological  sensor 
complement  of  the  Defense  Meteorological  Satellite  Program  (DMSP).  While  af- 
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fording  certain  advantages  characteristic  of  a  mature  technology  regarding  weight, 
power,  reliability,  and  cost,  passive  sensors  are  not  without  their  inherent  instru¬ 
mental  deficiencies.  Among  the.se  are  intrinsic  physical  limitations  on  attainable 
accuracy  and  vertical  resolution  and  the  inability  to  provide  operationally  useful 
measurements  of  some  desired  observables. 

In  order  to  effectively  fulfill  its  global  mission,  the  Department  of  Defense 
(DoD)  has  established  specific  observational  requirements  for  the  acquisition  of 
high  quality  geophysical  and  meteorological  data.  Among  the  stated  specifications 
for  each  desired  parameter  are  coverage  area,  horizontal  and  vertical  resolutions, 
mapping  accuracy,  and  the  range,  accuracy,  and  precision  of  the  measurement. 
Temporal  specifications  include  the  data  refresh  period  or  frequency  and  the  mea¬ 
surement’s  timeliness.  While  extant  satellite-borne  passive  sensors  partially  sat¬ 
isfy  some  of  these  needs,  remaining  gaps  in  operational  ob.servational  capabilities 
make  it  prudent  to  consider  alternative  remote  sensing  technologies.  Active  re¬ 
mote  sensor  systems  known  as  lidars  have  the  potential  to  bridge  some  of  these 
gaps,  particularly  those  due  to  the  inability  to  measure  a  given  parameter  using 
passive  means  or  deficiencies  in  accuracy  and  resolution.  Thus,  the  development 
of  space  qualified  laser-based  sensor  systems  for  actively  probing  the  atmosphere 
promises  to  significantly  enhance  the  state-of-the-art  of  satellite  remote  ssnsing. 

Numerical  weather  prediction  (NWP)  technology  and  practice  have  also 
evolved  significantly  since  their  operational  beginnings  in  the  1950s  and  1960s. 
Inferences  of  atmospheric  temperature  and  moisture,  clouds,  precipitation,  winds, 
and  surface  properties  are  currently  made  from  a  variety  of  passive  visible,  in¬ 
frared,  and  microwave  satellite  sensors.  For  global  NWP  it  is  primarily  the  re¬ 
trieved  temperature  profile  and  wind  data  which  are  used.  The  usefulness  of  other 
geophysical  parameters  for  NWP  is  not  well  established  and  is  largely  untested, 
although  a  great  many  retrieval  methods  have  been  proposed  or  developed  for  a 
variety  of  potentially  interesting  parameters.  Moisture  variables  -  that  is,  specific 
humidity,  clouds  and  precipitation  -  are  retrievable  and  are  potentially  very  useful, 
but  are  not  easily  assimilated  by  current  methods.  It  is  theoretically  possible  to 
retrieve  specific  htimidity  profiles  by  using  methods  analogous  to  those  used  to 
retrieve  temperatures.  However,  results  to  date  with  available  sensors  have  not 
been  wholly  satisfactory.  A  variety  of  other  meteorological  parameters  which 
may  be  used  for  NWP  are  retrievable.  These  include  winds  at  the  surface  from 
microwave  sensors  and  cloud  drift  winds  aloft,  as  well  as  other  surface  properties, 
such  as  soil  moisture,  albedo,  snow  cover,  temperature,  and  fluxes.  Positive  im¬ 
pact  of  current  satellite  observing  systems  on  global  NWP  has  been  demonstrated, 
especially  in  the  southern  hemisphere  where  conventional  data  sources  are  lim- 
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ited.  Increasing  the  accuracy  and  vertical  resolution  of  the  sensors  and  improving 
the  retrieval  and  assimilation  systems  should  lead  to  better  forecasts.  We  believe, 
nevertheless,  that  only  a  small  part  of  the  potential  usefulness  of  satellite  remote 
sensing  for  global  NWP  has  been  realized  and  that  lidar  .sensor  systems  may  play 
a  significant  role  in  this  evolution.  For  example,  previous  OSSEs  with  the  GL 
GDAS  (Hoffman  et  al.,  1990)  using  a  DWL  system  which  returns  complete  wind 
profiles  (WINDSA1')  have  shown  significant  positive  impact  upon  analyses  and 
forecasts  of  both  height  and  wind  relative  to  control  experiments.  Other  investi¬ 
gators  using  different  assimilation  .systems  have  also  found  WINDSAT  data  useful 
to  varying  degrees  (Arnold  ct  al.,  1985;  Atlas  ei  al.,  1985;  Dey  et  al.,  1985). 

Currently,  NASA  is  in  the  process  of  defining  the  instrument  and  platform 
characteristics  of  the  Laser  Atmospheric  Wind  Sounder  (LAWS)  as  part  of  the 
Earth  Observing  System  (EOS)  initiative.  As  pre.sently  defined  (Curran,  1989) 
the  LAWS  instrument  would  be  a  relatively  high  power  active  DWL  sensor  (10  J 
pulse)  operating  at  9  fim,  potentially  capable  of  returning  complete  profiles  of  hor¬ 
izontal  wind  from  the  top  of  the  troposphere  to  cloud  top  with  a  vertical  resolution 
of  approximately  1  km  and  a  horizontal  resolution  of  100  km.  Anticipated  accu¬ 
racies  are  on  the  order  of  1-3  ms“'  rms.  The  ba.seline  plan  has  LAWS  on  board  a 
polar  orbiting  satellite,  but  other  scenarios  are  under  study  (Atlas,  pers.  comm.). 

However,  there  are  legitimate  concerns  regarding  tlie  weight  and  power 
requirements  of  an  in.'itrumcnt  such  as  LAWS.  Secondly,  the  ability  of  the 
instrument  to  obtain  complete  vertical  profiles  of  wind  through  the  troposphere 
is  critically  dependent  upon  the  presence  of  sufficient  backscattering  aerosol. 
Data  from  the  Global  Backscatter  Experiment  (GLOBE)  should  help  resolve  this, 
but  the  answer  is  far  from  certain.  Finally,  the  size,  weight  and  power  of  the 
instrument  are  all  directly  proportional  to  its  ultimate  development  cost.  In  this 
regard,  it  has  recently  been  determined  that  LAWS  will  not  be  included  in  the 
instrument  package  slated  to  fly  aboard  the  Japanese  Pohu"  Orbiting  Platform 
(JPOP)  due  to  its  "high  risk"  nature. 

The  approach  adopted  for  this  study  is  to  assess  the  NWP  utility  of  a  down- 
scaled  (i.e.  lower  power  and  weight)  DWL  instrument  which  addresses  these, 
concerns  (note  that  the  primary  LAWS  objective  is  global  change  monitoring). 
The  scaled  back  instrument  would  be  smaller,  operate  at  lower  power  levels  (1  J 
pulse),  and  work  at  a  near-infrared  wavelength  (.see  Table  6).  These  characteristics 
would  lower  the  cost,  decrease  demands  upon  the  host  platform,  and,  relative  to  a 
comparable  infrared  lidar,  increase  the  back.scatter  signal.  Since  the  instrument  is 
operating  at  lower  jx)wer,  we  assutne  that  viable  return;:  will  only  occur  at  or  near 


"hard"  targets  where  backscattered  signal  is  siiflicient.  These  targets  are  cirrus 
cloud,  the  tops  of  other  clouds  and  marine  boundary  layer  aerosol. 

Finally,  we  note  that  range  gating  of  the  la.ser  pulse  it.self  can  allow  fairly 
accurate  estimation  of  the  target  height.  Although  not  utilized  in  the  experiments 
conducted  here  this  could  be  used  to  infer  cloud  top  height  which  might  also  be 
incorporated  into  a  data  assimilation  system. 

3.  TECHNICAL  APPROACH 

3.1  Lidar  Mca.suremcnt 

The  basic  principle  of  lidar  techniques  is  that  infonnation  about  the  physical 
state  of  the  atmosphere  or  hard  target  can  be  inferred  from  the  way  in  which 
the  laser  signal  is  backscattered  and  attenuated  by  atmospheric  constituents.  The 
theoretical  expression  relating  the  measured  signal  to  the  state  of  the  atmosphere 
is  the  lidar  equation: 

Er{u\  R)  -  Et{v)To{u*)T}{u')T^{u,u\  R)  ■ 

where  and  u'  are  the  transmitted  and  received  frequencies,  respectively.  Et 
is  the  transmitted  energy  and  Er  is  the  energy  received  from  a  range  cell  of 
thickness  AR  centered  at  a  distance  AR,  to  and  r/  are  the  system  parameters 
for  optical  transmission  and  filter  transmittance,  A  is  the  receiver  area,  is 
the  two-way  atmospheric  transmittance.  Pa  and  pm  are  aerosol  and  molecular 
backscatter  coefficients  respectively.  For  hard  target  measurements,  the  term  in 
brackets  containing  ilie  backscatter  quantities  is  replaced  by  surface  reflectivity  or 
albedo.  With  the  exception  of  winds,  information  on  the  state  of  the  atmosphere 
is  contained  in  the  tenns  /?,„  B,n  and  r^.  Wind  fields  can  be  derived  by 
measuring  the  Doppler  shift  of  the  back.scatter  signal  from  tracer  aero.sol  particles: 
A^'  =  i/'  -  u  =  2i>vlc  where  i/  is  the  component  of  wind  velocity  along  the 
sensor  line  of  sight. 

3.2  Doppler  Wind  Lidar  Sensors 

Aerosols  suspended  in  tlie  atmosphere  can  serve  as  wind  tracers  for  lidar 
measurements.  A  photon  back.scattered  by  an  aero.sol  panicle  moving  at  a  wind 
velocity  (z/)  in  the  line-of-sight  will  experience  a  Doppler  shift  in  frequency  of 
magnitude  {Au/u)  —  {2v/c),  where  c  is  the  speed  of  light.  Photons  .scattered 
by  particles  moving  toward  tlie  observer  in  tlie  linc-of-sight  will  experience 
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an  increase  in  frequency,  wliilc  tliosc  scattered  by  particles  receding  from  the 
observer  will  display  a  decrease  in  frequency.  'I’lie  Doppler  shift  caused  by 
aero.sol  backscatter  of  a  highly  stable  quasi-monochroinatic  laser  beam  could  be 
spectrally  analyzed  to  yield  i’’c  line-of-sight  component  of  wind  velocity.  In 
practice,  measuring  winds  by  thi.^  metluxl  is  very  difficult  becau.se  of  several  facts: 

•  The  Doppler  shift  is  extremely  small:  a  1  ms  '  wind  (line-of-sight)  result.'  in 
a  -3  X  lO  'nm  .shift  at  ()..'>  /on,  or  -6  x  10  -*'nm  at  10  //m. 

•  The  spacecraft  veUKity,  which  is  about  K  km  s  ',  also  contributes  to  the 
Doppler  shift.  The  contribution  of  this  velocity  along  the  line-of-sight  of  the 
measurement  must  be  detemiined  to  a  high  accuracy.  This  places  stringent 
requirements  on  instnimcnt  pointing  knowledge  and  spacecraft  attitude. 

•  The  aero.sol  back.scattcred  signal  is  often  quite  small  (particularly  in  the 
altitude  range  of  ."1-12  km  where  winds  are  very  importtmt)  and  at  some 
frc(|uencies  is  dwarfed  by  the  much  stronger  molecular  backscattcring  which 
acts  to  broaden  the  original  signal. 

There  are  two  major  techniques  using  coherent  and  incoherent  detection, 
to  determine  wind  velocity  in  this  way.  The  coherent  method  uses  heterodyne 
detection;  mixing  the  backscattered  .signal  with  a  local  oscillator  to  yield  a  beat 
frequency  proportional  to  the  Doppler  shift.  The  incoherent  technique  measures 
spectral  shifts  using  a  Fabry-Perot  interferometer  with  an  array-type  detector.  The 
relative  advantages  and  disadvantages  of  the  two  systems  are  discussed  by  Salvetti 
(1987),  and  Baker  and  Curran  (198.5).  llicse  tlircc  works  also  summarize  potential 
implementations  of  DWLs.  It  is  now  considered  technically  feasible  to  measure 
wind  profiles  from  space  by  measuring  the  Doppler  shift  of  a  transmitted  laser 
pulse  (Salvetti,  1987;  Curran  and  others,  1988).  With  a  strong  enough  signal  the 
reflected  pulse  may  be  range  gated  to  yield  vertical  resolution  of  1  km  or  better. 
The  strength  of  the  reflected  signal  depends  principally  on  the  energy  tran.smitted 
and  the  reflectivity  of  the  atmospheric  volume  being  sampled.  Since  the  laser 
may  be  ftKuscd  to  a  ver>’  fine  solid  angle,  horizontal  resolution  as  fine  as  desired 
(down  to  .scales  of  meters)  may  be  obtained.  Further,  the  lidar  measurements  are 
a  direct  measure  of  the  line  of  sight  (or  radial)  vel(x:ity.  Thus  if  the  reflected 
signal  strengths  are  suflicient  any  reasonable  desired  accuracy  might  be  obtained 
by  this  technique.  Two  measurements  of  the  same  atmospheric  volume  from 
different  viewing  angles,  along  with  the  assumption  that  vertical  velocities  are 
negligible,  arc  required  to  infer  the  u  and  v  wind  components.  The  simplest 
method  of  accomplishing  this  is  to  use  a  conical  .scan  pattern  (cf.  Fig.  19  of 
Curran  ct  al.  (1988)).  The  global  distribution  of  aero.sol  is  not  well  known.  Since 
tl.e  aero.sol  concentration  directly  effects  the  atmospheric  reflectivity  to  the  lidar 


signal,  the  relationship  of  tlic  DWL  errors  to  the  transmitted  energy  cannot  be 
reliably  predicted.  I  lowevcr,  in  the  experiments  conducted  here,  we  obtained  an 
estimate  of  measurement  error  by  assuming  .several  aerosol  and  cloud  drop  size 
distribution  scenarios.  In  particular,  we  assume  sufficient  aerosol  particles  in  the 
marine  boundary  layer  to  enable  a  measurement  to  be  made. 

3.3  NWP  Related  Lidar  In.strunient  Characteristics 

From  the  point  of  view  of  numerical  weather  prediction  (NWP),  the  most 
important  characteristics  of  any  jiroposcd  remote  sensing  system  are  its  geograph¬ 
ical  coverage,  horizontal  and  vertical  resolution  and  its  error  cinuacteristics.  In 
a  simulation  study  these  characteristics  must  be  properly  accounted  for.  'ITiesc 
considerations  lead  immediately  to  a  number  of  issues  which  bear  on  simulation 
experiments. 

1.  Lidar  coverage  depends  on  cloudiness.  Since  the  lidar  pulse  penetrates  the 
cloud  only  weakly  and  since  multiple  scattering  from  the  topmost  part  of  the 
cloud  contaminates  any  signal  received  from  below,  it  is  impossible  to  retrieve 
wind  profiles  below  clouds.  In  the  present  experiments,  there  will  be  no  lidar 
data  at  points  below  large  scale  cloudiness.  This  is  perhaps  pessimistic,  since 
on-board  shot  management  might  find  holes  in  some  of  the  cloud  fonnations. 

2.  For  NWP,  it  is  not  just  accuracy  of  the  measurement  which  is  imponant, 
the  measurement  must  be  repre.sentative  as  well.  NWP  is  really  concerned 
with  the  spatially  and  temporally  smoothed  behavior  of  the  atmosphere.  That 
is  variations  on  the  .scale  of  nteters  and  .seconds,  in  fact  on  the  scale  of 
kilometers  and  minutes,  are  ge-  crally  considered  to  be  averaged  over  and 
are  parameterized  witliin  the  model.  Consequently,  that  part  of  the  measured 
signal  attributable  to  these  .scales  is  considered  to  be  noise  from  the  NWP 
point  of  view.  This  source  of  error  can  in  some  cases  be  predominant.  The 
representativeness  issue  is  of  some  importance  to  lidar  measurements.  Lidars 
naturally  average  in  the  vertical.  This  is  gootl  for  NWP,  however,  this  average 
is  weighted  by  the  venical  profile  of  aerosol  concentration.  Our  concern  lies 
in  the  horizontal  shot  pattern.  .Since  the  atmospheric  volume  sampled  may  be 
only  meters  across  and  there  is  no  tcmitoral  averaging,  lidar  measurements 
will  contain  large  errors  of  representativeness.  Averaging  many  shots  would 
oscicome  this  error  source,  but  shot  patterns  which  have  been  propo.scd  have 
a  typical  separation  between  neighboring  shots  of  100  km. 

3.  Spatially  correlated  errors  are  difficult  for  an  analysis  .scheme  to  remove, 
because  the  ckita  tend  to  corroborate  each  other.  Real  data  tend  to  itave 
correlated  errors.  Fven  for  radiosondes,  significant  vertical  error  correlations 
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are  present.  DWL  data  errors  might  be  strongly  correlated,  for  example,  if 

they  are  due  to  an  error  in  the  estimate  of  spacecraft  vekKity. 

4.  OSSE  DESIGN 

llie  OSSli  which  was  conducted  is  compo.sed  of  four  components:  (i)  the 
nature  run,  a  simulation  of  the  atmosphere,  (ii)  a  simulation  of  the  observing 
system,  (iii)  a  quasi-operational  global  NWP  4D  assimilation  and  forecast  system, 
and  (iv)  a  comparison  metluKlology.  Many  aspects  of  these  components  have 
been  described  at  length  in  previous  reports  and  publications  (Hoffman  ct  al., 
1990;  Louis  et  at.,  1989)  and  arc  only  briclly  di.scussed  here. 

4.1  The  Nature  Run 

The  nature  run  is  the  data  .set  from  which  simulated  data  are  extracted.  Ideally 
the  nature  run  would  be  a  detailed  history  of  the  real  atmosphere.  Unfortunately, 
any  series  of  analy.ses  of  the  atmosphere  contains  substantial  errors  due  to  ob¬ 
servational  errors,  analysis  methodology  errors,  forecast  model,  and  initialization 
errors.  Because  of  this,  O.SSLs  reipiire  a  simulated  nature  run.  'Riis  model  should 
be  different  from  the  model  used  in  the  4D  assimilation  to  avoid  the  identical 
twin  problem  in  which  the  simulated  analy.ses  are  unrealistically  similar  to  tlie 
nature  run.  This  model  must  also  have  high  resolution,  sophisticated  physics,  and 
a  reasonable  climtitology. 

In  this  study,  the  nature  run  is  a  20  day  ECMWF  forecast  beginning  at 
(X)  UTC  10  November  1979  from  the  Global  Weather  E.xperimcnt  (GWE)  Illb 
analysis  produced  at  the  European  Center  for  Medium  Range  Weather  Forecasting 
(ECMWF)  (Bengts.son  et  at.,  1982).  The  mtxlel  used  in  the  nature  run  forecast 
was  a  version  of  the  15  layer,  1.875  degree  gridpoint  model  (Hollingsworth  et  at., 
1980).  This  model  included  fairly  complete  physics  (Ticdtke  et  at.,  1979)  with  a 
diurnal  cycle.  We  used  the  archived  forecast  fields  which  were  available  every 
12  hours,  intcr|iolated  to  2.5  degree  resolution.  Since  the  nature  run  is  relatively 
smooth,  this  slightly  degraded  resolution  is  a  satisfactory  representation. 

4.2  The  Sinuilated  OI)scrving  Systems 

'Hie  simulated  observing  sy.stems  should  be  as  sophisticated  as  possible  with 
reasonable  computation  constraints.  Characteristics  of  all  observations  and  their 
errors  must  be  treated.  In  particular,  with  respect  to  the  lidar  data,  the  clouds  and 
aerosols  predicted  by  the  nature  mtxlel  should  be  used  in  determining  the  spatial 
distribution  of  the  data  and  its  overall  accuracy.  Note  that  the  most  interesting 
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meteorology  is  often  associated  with  clouds.  The  horizontal  and  vertical  sampling 
characteristics  of  the  wind  sensor  must  be  properly  modeled.  Here  we  briefly 
de.scribe  the  other  simulated  data  which  were  used  in  the  OSSEs.  Details  of  the 
lidar  simulation  methodology  are  contained  in  the  next  section. 

With  the  exception  of  the  lidar  ob.servations,  all  data  were  simulated  at  the 
National  Meteorological  Center  (NMC)  by  replacing  nearly  all  valid  GWE  level 
Ilb  observations  obtained  during  the  special  ob.serving  periods  with  simulated 
observations  at  'he  same  IcKations.  The  ob.serving  systems  include  conventional 
observations  such  as  raobs  and  aireps,  as  well  as  satellite  cloud  drift  winds 
tCDW),  and  satellite  temperature  soundings  (TOYS).  The  simulated  standard  data 
were  created  by  interpolating  nature  data  to  the  GWE  obsers'ation  and  adding  a 
random  observing  error.  Tliis  procedure  yields  quite  reasonable  data  coverage 
and  quality  control  in  the  simulated  data.  However,  certain  discrepancies  are 
possible:  for  example,  cloud  drift  winds  may  be  present  where  there  are  no  clouds 
in  the  nature  run.  With  regard  to  the  CDWs,  these  simulated  ob.servations  were 
treated  differently  from  the  CDWs  used  in  previous  OSSEs.  The  errors  of  these 
observations  were  modified  to  be  more  representative  of  real  CDW  measurements, 
as  de.scribed  below. 

4.2.1  Filtering  of  CDW  Errors 

The  CDWs  simulated  by  Dey  ct  al.  (1985)  have  random  error  only,  whereas 
we  know  a  sizeable  fraction  of  the  enor  is  due  to  height  assignment  errors 
(McPherson,  1984).  Eurther  sources  of  error  arc  the  fact  that  cloud  ba.se,  not  the 
cloud  top,  may  be  the  bc.st  level  for  assigning  the  winds  and  the  deviations  from 
unit  emissivity  for  the  cloud.  These  errors,  especially  the  height  errors,  tend  to 
be  very  well  correlated  (at  least  at  a  particular  synoptic  time)  for  a  particular  data 
priKluccr.  The  cau.se  of  the  height  assignment  errtirs  is  fundamental.  Generally 
two  approaches  to  height  assignment  have  been  used: 

1.  The  climatological  ap|uoach.  For  example,  assigning  all  low  level  wind  to 
9(X)  mb  or  assigning  all  high  level  wind  to  the  climatological  tropopau.se. 

2.  IR  radiance  approach.  In  this  technique  the  ob.scrved  IR  cloud  radiance  is  used 
to  deduce  a  cloud  toji  temjierature,  which  is  then  matched  to  a  temperature 
profile  specilicd  from  a  recent  analysis  or  short  term  forecast. 

The  CDWs  simulated  by  Dey  ct  al.  essentially  u.sed  an  error  free  temperature 
profile  for  height  assignment.  We  could  reassign  the  pressures  by  adding  an  ad 
hoc  error  to  the  CDW  pressures  directly,  or  by  assuming  an  ad  lux:  temperature 
error  and  reassigning  the  pressure  by  comparing  the  “observed”  temperature  to 
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the  nature  run.  The.se  errors  would  have  to  be  fairly  large.  These  approaches 
would  indirectly  create  horizontally  correlated  errors. 

We  have  cho.sen  to  more  directly  induce  correlated  errors  in  the  CDW  data 
by  spatially  filtering  the  errors  in  the  simulated  data.  We  divided  the  data  into 
high  and  low  winds  and  filtered  the  u  and  v  wind  components  for  each  group 
.separately.  As  seen  in  figure  1  the  distribution  of  oh.serving  pressures  has  two 
major  modes  and  we  chose  6()0  mb  as  the  dividing  pressure.  The  data  are  for  12 
UTC,  21  November  1979;  there  were  a  total  of  1935  CDWs.  The  filter  we  use  is 
a  weighted  average  followed  by  :in  amplification  step.  The  averaging  reduces  the 
amplitudes  of  the  original  errors,  while  inducing  correlations.  The  amplification 
step  multiplies  all  weighted  average  errors  by  a  constant  factor  to  restore  the 
original  amplitude  of  the  error  variance,  lliis  constant  is  equal  to  the  variance 
of  the  original  errors  divided  by  the  variance  of  the  filtered  errors.  The  weighted 
average  operator  is  equivalent  to  the  first  pass  of  a  Barnes  filter  (Koch  et  al., 
1983).  Hie  weights,  h-,  are  given  in  terms  of  the  distances,  d,  between  the  point 
to  be  filtered  and  the  nearby  points,  according  to  w  =  cxp{-(dldc,  Here  do  is  the 
scale  of  the  filter  and  all  points  within  2xdo  are  included  in  the  weighted  average. 
After  filtering  a  random  error  field,  correlations  drop  to  approximately  0.5  at  a 
separation  distance  of  do.  We  chose  do  to  be  0. 1  earth  radii  or  approximately 
637  km. 

'file  results  of  the  filter  on  the  horizontal  correlations  of  the  CDW  errors  is 
shown  in  figures  2  and  3  for  the  low  and  high  wind  groups.  For  reference  the 
correlation  of  the  nature  nin  evaluated  at  the  .same  data  locations  is  presented  in 
the  upper  row  of  plots.  Ccrrclations  of  the  original  errors  are  displayed  in  the 
center  row  and  correlations  of  the  filtered  errors  are  displayed  in  the  lower  row 
of  plots.  As  indicated  the  results  for  the  u  and  v  wind  components  are  displayed 
in  the  left  and  right  columns  respectively.  Note  that  the  errors  are  treated  by  the 
filter  as  if  they  were  all  at  a  single  level  (.separately  for  high  and  low  winds). 
On  the  other  hand,  the  nature  data  used  for  the  correlations  are  at  a  variety  of 
pressure  levels,  'fhe  .synoptic  scale  is  clearly  evident  in  the  v  correlations  from 
the  nature  run.  This  is  not  obvious  in  the  u  correlations  because  of  the  strong 
zonal  component,  especially  at  higher  levels. 
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Figure  2.  Horizontal  correlations  fur  low  level  (1(HX)-(>(K)  mb)  CDWs.  Shown  are  correlations 
for  tlie  backgrminil  field  evaluated  at  the  C’DW  locations  (top),  original  CDW  cnors 
(middle)  .iiid  filtered  CDW  ctrors  (Nittoin)  lor  u  (lell)  and  v  (right)  wind  components. 
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Figure  3.  liori/onial  coiicl.iiitins  for  high  level  (WKI-MK)  mb)  ClW.s.  Same  as  Fig.  2. 
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4.3  The  41)  Assimilation  System 

Each  simulated  data  assimilation  experiment  consists  of  one  assimilation  run 
from  00  UTC  18  November  through  00  UTC  25  Novem.ber  1979  and  three  96- 
hour  forecasts  from  00  UTC  21,  23,  and  25  November.  Each  assimilation  run 
consists  of  a  series  of  assimilation  cycles,  and  each  cycle  in  turn  is  made  up  of  a 
6-hour  forecast  that  serves  as  the  background  or  first  guess  field  for  the  analysis, 
an  optimal  interpolation  analysis  that  combines  the  first  guess  fields  with  the 
observations,  and  a  nonlinear  normal  mode  initialization  of  the  analysis.  The 
initialized  analysis  is  the  starting  point  for  the  next  6-hour  forecast,  which  again 
serves  as  the  first  guess  field  for  the  sub.sequent  analysis. 

The  AFGL  Statistical  Analysis  Program  (ASAP)  (Norquist,  1986;  Norquist, 
1988)  was  developed  from  the  NMC  multivariate  optimal  interpolation  (01)  proce¬ 
dure  as  described  by  Bergman  (1979)  and  by  McPherson  et  cl.  (1979).  The  01  is 
a  multivariate  analysis  of  height  and  wind  components  and  a  univariate  analysis  of 
relative  humidity.  The  corrections  for  an  analysis  grid  point  arc  weighted  sums  of 
surrounding  ob.servation-minus  first-guess  residuals.  Equatorward  of  70  degrees, 
these  weights  are  calculated  as  described  by  Dey  and  Morone  (1985).  In  the  polar 
regions,  the  Bergman  (1979)  formulation  is  used.  The  analysis  is  performed  in  the 
sigma  coordinates  of  the  model  on  a  Gaussian  grid  of  62  x  61  latitude-longitude 
points.  Recent  changes  to  the  01  and  documentation  are  discussed  below. 

The  GL  normal  mode  initialization  (NMl)  is  based  on  the  NMC  NMI  (Ballish, 
1980).  The  GL  global  specU'al  model  is  based  on  the  NMC  GSM  designed 
by  Sela  (1980).  For  the  version  used  here,  the  physics  routines,  taken  almost 
intact  from  NMC  (circa  1983),  include  surface  exchanges  of  heat,  moisture, 
and  momentum  over  the  ocean;  surface  exchanges  of  momentum  over  land;  a 
dry  adiabatic  adjustment;  large  scale  supersaturation  precipitation;  and  a  Kuo 
convection  scheme.  The  hydrodynamics;  i.e.,  the  adiabatic,  inviscid  dynamics 
including  vertical  and  I'  )rizontal  advection,  time  stepping,  and  transformations 
between  spectral  and  physical  space;  were  redesigned,  as  documented  by  Brenner 
etal.  (i982;1984). 

Tlie  spectral  resolution  of  the  forecast  model  itself  is  defined  by  a  rh  mboidal 
truncation  at  wave  number  30.  The  Gaussian  grid  of  the  forecast  model  contains 
76  X  96  latitude-longitude  points.  There  are  12  vertical  layers,  the  top  5  of  which 
have  no  moisture. 


4.4  Comparison  Methodology 

To  determine  possible  impact  of  the  lidar  observing  system,  output  from  the 
GDAS  is  examined  for  cases  run  with  (LIDAR)  and  without  (CONTROL)  the 
inclusion  of  the  DWL  data.  Table  1  summarizes  data  usage  for  the  CONTROL 
and  LIDAR  experiments.  Subjective  evaluations  are  made  by  looking  at  the 
standard  meteorological  variables  (e.g.  horizontal  wind,  gcopotential  height,  and 
relative  humidity)  on  synoptic  charts.  More  objective  comparisons  arc  also  made 
using  calculations  of  global  and  regional  error  statistics. 


T:ibIo  1.  OSSE  l)al.T  Usage. 


Conventional 

CDW 

TOYS 

DWL 

CONTROL 

X 

X 

X 

LIDAR 

X 

X 

X 

X 

Tlie  CONTROL  experiment  represents  a  new  baseline  case  which  was  run  to 
measure  the  DWL  impact.  It  differs  from  the  original  baseline  experiment  known 
as  STATSAT  used  in  earlier  OSSEs,  (Hoffman  et  al.,  1990)  in  that  it  features 
more  realistic  CDWs  (see  above)  and  a  modified  01. 

Essentially,  changes  to  the  analysis  system  were  made  in  three  main  areas: 

•  Vurhible  resolution 

A  series  of  changes  were  made  to  the  analysis  ctxles  (ASAP  and  MASAP) 
and  their  associated  pre-  and  postprocessors  to  allow  for  variable  resolution. 
The  changes  to  the  analysis  codes  were  made  as  a  series  of  small  changes,  so 
that  the  results  of  the  initial  code  with  small  errors  could  be  reproduced,  the 
source  of  which  could  be  clearly  identified.  Significant  changes  have  been 
reflected  in  the  new  documentation  for  ASAP  (August  1989).  Note  that  all 
common  blocks  are  now  "included"  and  arc  the  same  whenever  they  are  used. 
The  constants  for  LOWTMP  (the  1-lattery  algorithm)  are  calculated  once  at 
the  start  of  ASAP,  since  these  constants  depend  on  the  vertical  structure.  The 
code  is  now  .strictly  sttindard  1-ORTRAN.  The  largest  differences  during  this 
series  of  mcxiili  cat  ions  :irc  due  to  the  calculation  of  the  LOWTMP  constants  to 
full  machine  precision  instead  of  specifying  them  on  a  DATA  statement  to  9  or 
10  digits.  Modifications  for  MASAP  paralleled  tho.se  for  ASAP.  In  particular 
all  commons  are  now  "included".  Previously  the  finst  guess  error  growth  rates 
on  the  nuKiel  sigma  structure  were  hardwired  in  the  cixle.  Now  the.se  rates  are 
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stored  at  mandatory  levels  and  interpolated  to  the  sigma  structure.  This  was 
the  one  change  which  produced  errors  detectable  at  the  0(10'^  g  kg  ’)  level. 

'Hie  pre-  and  postprocessors  (PPP.s)  were  completely  redesigned,  in  a  modular 
fashion.  Since  the  GSM  PPPs  already  exist  in  variable  resolution  fonnat,  we 
only  implemented  our  design  for  the  PPPs  for  the  01. 

•  Surface  pre.ssure  analy.si.s 

Tlie  surface  pressure  update  was  reformulated.  The  previous  update  was 
described  by  Norquist  (1986,  pp.  122-123)  and  involved  extrapolating  the 
heights  at  the  top  of  the  lowest  three  layers  for  both  the  background  and 
analysis  to  the  a  -  \.  The  new  formulation  u.ses  the  lowest  and  underground 
layer  height  increments  to  e.stiniate  the  height  increment  at  cr  =  1,  which 
is  then  converted  hydrostatically  to  a  pre.ssure  increment.  This  approach  is 
similar  to  that  used  by  the  ECMWF. 

•  Use  of  RAOB  hcigii(.s 

The  radiosonde  heights  interpolated  to  the  sigma  levels  are  now  used  directly 
as  observations.  Previously,  these  data  were  u.scd  only  for  thicknesses;  they 
were  differenced  hydrostatically  to  obtain  layer  temperatures  which  were  then 
anchored  to  the  first  guess  heights  (usually  at  the  top  of  layer  1)  or  to  the 
mtxlel  topography.  The  latter  case  would  only  hold  if  the  raob  profile  started 
below  the  model  topography,  while  in  the  fonner  case  the  first  height  datum 
would  be  at  the  top  of  the  layer  above  the  anchoring  level  (normally  layer  2). 
Consequently,  this  change  increa.sed  the  number  of  raobs  at  the  lowest  levels 
as  indicated  in  table  2  below. 


Tabic  2.  Number  of  Kadio.sondc  Ob'arrvation':  at  the  Top  of  the 
First  Four  Sipina  l..•ly^TS  for  1200  UTC  21  November  1979. 


a  Ixvel 

1 

2 

3 

4 

.ST.ATSAT 

2,^0 

4^5 

671 

679 

CON'l'ROL 

439 

673 

690 

679 

The  use  of  ab.solute  heights  is  probably  more  imptfrtant  than  the  increase  in 
the  number  of  low  level  observations. 
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These  changes  were  incorjx)rated  incrementally  to  isolate  the  relative  impacts 
of  each,  and  to  detenu ine  the  effect  of  machine  word  length  on  the  analyses 
(testing  was  done  both  on  the  GL  Cyber  and  the  AER  Alliant).  All  the  above 
mtxlifications  have  been  described  at  length  in  the  contract  quarterly  status  reports 
(no.  2,  no.  3,  and  no.  4;  GL  Contract  I'19628  -89  -C-(K)44). 

In  Figure  4  rins  height  errors  are  shown  for  both  CONTROL  and  Sl'ATSAT 
assimilation  runs.  These  indicate  a  substantial  reduction  in  analysis  error  of 
approximately  10  m  for  CONTROL.. 

5.  LIDAR  DATA  SIMULAFION  METHODOLOGY 

The  methodology  used  to  simulate  the  lidar  data  includes  a  sampling  strategy 
(i.e.  where  to  locate  potential  measurements),  a  cloud  and  aerosol  diagnosis  (since 
these  determine  the  vertical  and  horizontal  distribution  of  viable  lidar  returns),  and 
measurement  errors  (which  arc  a  function  of  instrumental  characteristics,  scanning 
geometry,  and  assumed  cloud/aerosol  backscatter  characteristics). 

5.1  Data  Sampling 

As  noted  above,  an  extended  20  day  forecast  of  the  ECMWF  grid  point 
model  serves  as  our  ground  truth  for  verification  and  cloud  diagnosis.  This  is  the 
same  nature  run  used  in  previous  OS.SEs  (Hoffman  eiai,  1990).  The  data  (wind, 
temperature,  and  relative  humidity)  on  a  2..*)  degree  re.solution  grid  are  intcrjiolated 
in  space  and  time  to  the  horizontal  legation  of  the  anticipated  lidar  measurement. 

The  actual  locations  of  the  lidtu-  measurements  have  been  determined  using 
the  orbital  parameters  of  the  DMSP  polar  orbiting  platform  and  the  instrumental 
scan  parameters  of  the  cross-tracking  Special  Sensor  Microwave  temperature 
(SSM/T)  sounder  (Isaacs  a  oL,  198.i).  In  practice  the  lidtu'  instrument  will  be 
conically  scanning  with  each  of  the.se  prescribed  kx'ations  the  result  of  intersecting 
fore-  and  aft-looking  shots  (line  of  sight  measurements)  binned  in  a  grid  like 
that  of  the  SSM/T  instrument.  Thus,  the  important  geometric  features  of  the 
conically  .scanning  lidar  may  be  deduced  from  the  observation  locations.  The 
scanning  geometry,  in  turn.  alTeels  the  nature  of  tlie  measurement  errors,  which 
are  discussed  below. 
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5.2  Cloud  and  Aerosol  Diagnosis 

In  our  simulations  three  types  of  lidar  returns  are  allowed  for;  (1)  liquid  water 
cloud,  (2)  cirrus  cloud,  and  (3)  marine  boundary  layer  aerosol. 


5.2.1  Liquid  Water  Cloud 


Liquid  water  clouds  :ue  diagno.scd  using  the  Gcleyn  .scheme.  In  this  formu¬ 
lation  the  cloud  fraction  at  each  nature  run  moist  level  is  parameterized  in  terms 
of  the  actual  relative  humidity  and  a  predefined  critical  relative  humidity  at  that 
level  as: 


CO(k)  = 


'  RHik)  -  RHc{k) 

100.  -  RHcik) 


2 


(5.1) 


wh(.»e  RH(k)  is  the  relative  humidity  at  level  k,  RHc(k)  is  the  critical  relative 
humidity  at  level  k,  and  CC(k)  is  the  computed  cloud  fraction  at  level  k.  RHc  is 
itself  a  function  of  pressure  and  we  have  u.scd  the  same  values  that  were  used  in 
the  cloud  parameterization  scheme  during  the  forward  integration  of  the  nature 
run.  Values  of  RHc  at  the  six  pressure  levels  from  1(XX)  mb  to  3(X)  mb  are:  100.0, 
59.0,  50.0,  60.3,  72.5.  Note  that  this  diagnosed  cloudiness  is  only  the  large  scale, 
non-convective  amount. 


5.2.2  Cirrus  Cloud 


Cirrus  clouds  aie  diagnosed  from  both  the  relative  humidity  (RH)  and  tem¬ 
perature  (T)  at  the  ob.scr\ing  location.  We  have  adopted  a  simple  thresholding 
approach  such  that  a  cirrus/ice  cloud  is  present  whenever  both  the  T  and  RH 
fall  below  and  above  their  critical  values,  respectively.  The  critical  T  has  been 
fixed  throughout  at  253  K,  although  a  dependence  on  latitude  and  level  mtght  be 
considered.  A  latitude-dependent  critical  RH  was  used  to  simulate  this  feature 
of  published  cirrus  climatologies.  Values  of  critical  RH  were  tuned  to  approxi¬ 
mate  the  fall  time  peritxl  zonal  climatology  derived  from  Stratospheric  Aerosol 
and  Gas  Experiment  (SAGE)  data  shown  by  Woodbury  and  McCormick  (1986). 
These  values  are  slntwii  below. 
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Tabic  3.  Critical  RH  Values  for  Cirrus  Cloud. 


Latitude  (deg.) 

RHc(%) 

-20  -»  +20 

60. 

±  20  -»  ±  40 

80. 

.+  40  -r  +  70 

70. 

±  70  ±  90 

85. 

Comparison  with  the  SAGK-derived  climatology  is  shown  in  figure  5.  These 
results  correspond  to  one  six  hour  interval  during  the  7  day  assimilation  period, 
but  are  generally  representative  of  other  times  during  the  period  as  well. 

Since  the  nature  run  data  only  contains  moisture  information  up  to  300  mb 
we  need  to  allow  for  the  (Kcurrence  of  cirrus  at  higher  levels.  To  do  this  we 
assume  that  any  time  cirrus  is  diagnosed  at  300  mb  it  extends  all  the  way  up  to 
the  tropopau.se  level.  The  tropopause  level  is  a  fixed  function  of  latitude  ranging 
from  100  mb  in  the  tropics  to  3(X)  mb  in  the  polar  regions  as  indicated  in  Table 
4  below.  No  attempt  is  made  to  deduce  a  cirrus  cloud  fraction  or  cloud  opacity 
from  the  nature  run. 


Table  4.  Laiiiudc  Dependent  Tropopause  Pressure  used  in  Cirrus  Diagnosis. 


Latitude  (degrees) 

P  U'opopause(n^fi) 

±  70  -»  ±  90 

300. 

±50  ->±  70 

250. 

±  30  •  »  ±  .50 

200. 

±  10  ■  >  ±  .30 

1.50. 

0  ±  90 

KM). 

Cloud  frequency 


latitude 


Figure  5.  Qimpjrison  ol  cinus  occurrence  dclcnnincd  from  SAGE  dala  vs. 
dial  obtained  in  currenl  O.SSE.  SAGil-  resiills  valid  for  Sepienitier,  Oclobor, 
November  lime  period.  AEK  roMillN  valid  for  18  Novemlx'r  1979.  06  LITC. 
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5.2.3  Marine  Aerosol 

We  have  assumed  that  the  marine  boundary  layer  will  always  contain  sufficient 
aerosol  for  a  lidar  measurement.  Therefore,  any  observing  location  over  the 
unfrozen  (Keans  (sea  ice  extent  is  .set  at  November  climatology)  contains  marine 
aerosol  in  the  boundary  layer,  which  is  assumed  to  be  of  fixed  depth  with  a  top 
at  850  mb. 

5.2.4  Measurement  Decision  Algorithm 

The  prtKess  whereby  the  particular  geophysical  profile  at  an  observing  IcKa- 
tion  is  analyzed  for  the  levels  of  viable  lidar  returns  is  essentially  a  decision  tree. 

After  a  profile  is  analyzed  for  the  liquid  water  cloud  fraction  (if  any)  at  each 
level,  the  profile  is  integrated  downward  from  the  top  of  the  atmosphere  (assuming 
random  overlap)  to  obtain  integrated  cloud  amounts  at  each  level. 

Now  the  profile  may  be  inspected  level  by  level  from  the  top  of  the  atmosphere 
downward  and  a  return  type  assigned  to  each  level.  Each  level  is  checked  in  turn 
for  liquid  water  cloud,  cirrus,  and  then  marine  aerosol.  For  the  case  of  cirrus  cloud, 
its  presence  is  sufficient  to  permit  a  measurement  since  no  information  on  cloud 
fraction  or  opacity  is  known.  Additionally,  cirrus  is  assumed  to  be  sufficiently 
non-opaque  so  that  its  pre.sencc  does  not  preclude  returns  from  lower  level  water 
clouds  or  aerosol.  For  liquid  water  clouds  and  marine  aerosols  a  measurement 
(Kcurs  if  the  probability  that  the  cloud  level  or  aerosol  layer  is  seen  exceeds  a 
critical  value.  This  probability  is  just  the  fraction  of  the  layer  or  level  seen  from 
above  and  is  computed  directly  from  the  integrated  cloudiness  profile.  These 
critical  value  thiesholds  are  important  since  they  affect  the  efficiency  with  which 
partially  visible  levels  yield  wind  measurements  and  arc  inversely  proportional 
to  the  assumed  strength  of  the  backscattered  signal.  Therefore,  during  the  data 
simulation  critical  values  for  water  cloud  and  aero.sol  probabilities  were  set  at 
0.3  and  0.5,  respectively,  since  cloud  water  droplets  have  stronger  scattering 
characteristics  at  near  infrared  wavelengths.  Table  5  indicates  the  lidar  return  type 
totals  by  level  for  one  six  hour  time  interval  during  the  assimilation  period  in  which 
4044  profiles  were  processed.  In  the  table  the  tenns  warm  and  cold  cirrus  refer 
to  cirrus  (Kcurring  at  temperatures  above  and  below  -40  degrees  C,  respectively. 
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According  to  Heymsfield  and  Platt  (1984)  this  temperature  corresponds  to  a  sharp 
change  in  observed  ice  crystal  size  distribution  and  lidar  backscatter  regimes.  At 
the  lower  temperatures  the  ice  particle  concentrations  are  seen  to  decrease  along 
with  lidar  backscatter. 


Table  5.  I.idar  Kelurn  Totals  by  Typo  and  Ixjvel  for  6  Hour  Period. 


Pressure 

(mb) 

Total 

Returns 

Water 

Cloud 

Wann 

Cirrus 

Cold 

Cirrus 

Marine 

Aerosol 

50. 

0 

0 

0 

0 

0 

70. 

0 

0 

0 

0 

0 

100. 

192 

0 

192 

0 

0 

l.iO. 

432 

0 

432 

0 

0 

200. 

457 

0 

457 

0 

0 

250. 

457 

0 

457 

0 

0 

3(K). 

1287 

0 

4.57 

8.30 

0 

4(K). 

1 101 

14 

.564 

613 

0 

500. 

1226 

21 1 

794 

221 

0 

7(M). 

1229 

710 

519 

0 

0 

850. 

735 

735 

0 

0 

0 

MBL 

1401 

0 

0 

0 

1401 

Figure  6  shows  tlie  locaiion  of  lidar  measurements  at  8.'i()  mb  for  the  6  hour 
period  centered  at  12  IJTC  on  19  November.  The  ob.scrvaiions  shown  correspond 
to  areas  <rf  low  level  cloudines.s.  Note  again  that  the  distribution  of  observations 
i.^  far  from  conii)leie. 


5.3  Measurement  Error 

Once  the  return  types  have  been  set  for  each  vertical  level  the  data  are  created 
by  adding  an  error  to  tlie  "exact"  u  and  v  components  at  each  level  found  to  have 
a  return.  The  "exact"  values  are  obtained  directly  from  the  nature  run  profile 
with  one  exception:  the  wind  in  the  marine  boundary  layer  is  considered  to  be  an 
average  of  the  wind  components  at  1000  mb  and  850  mb. 

5.3,1  Cloud  Optical  Properties 

Cloud  optical  propeiiics  are  an  iiuportanl  issue  because  they  determine  the 
strength  of  the  backscatier  signal  received  by  the  sen.sor  and  ultimately  tlic 
accuracy  of  the  measured  wind.  Spec  fically,  we  desire  an  estimate  of  tlie  cloud 
backscatter  for  various  cloud  types.  For  liquid  water  clouds  and  boundary  layer 
aerosol  we  adopted  well  known  models  for  stratus  and  stratocumulus  drop  size 
distributions  and  marine  boundary  layer  aerosol  distributions.  The  estimation  of 
cirrus  optical  properties  including  extinction  and  the  desired  lidar  backscatter  cross 
sections  was  based  on  the  parameteiization  of  the  cirrus  particle  size  spectrum 
by  Heymsfield  and  Platt  (1984).  This  parameterization  provides  an  empirically 
based  relationship  between  the  particle  size  spectrum  of  ice  clouds  and  the  ambient 
temperature  and  ice  water  content.  The  ice  crystal  size  distribution  over  the  size 
range  D  >20  fim  was  represented  by  two  power  law  equations  of  the  form: 

n{D)  =  IWC  •  Ai  •  for  t=  / ,  2  (5.2) 

where  n  has  the  dimensions  of  m  D  is  the  maximum  crystal  dimension 

in  /xm,  and  IWC  is  the  ice  water  content  in  gm  -'.  Hie  subscript  refers  to  the 
break  point  in  the  distribution  given  by  a  crystal  maximum  dimension  Do.  The 
constants  .i,  and  Bi  are  related  to  the  ice  concentrations  at  fixed  sizes  and  the 
size  distribution  slopes,  respectively.  Values  of  the  con.stants  and  the  ice  water 
content  as  a  function  of  ambient  temperature  in  the  range  of  -20  to  -60  degrees 
C  are  given  in  Table  2  of  Heymsfield  and  Platt  (1984).  As  noted  above,  one 
result  of  their  analysis  is  that  an  apparent  change  in  ice  particle  size  distribution 
occurs  between  -35  and  -40  degrees  C.  Tltis  results  in  a  different  class  of  size 
distributions  corresitonding  to  "warm"  and  "cold"  cirrus,  respectively.  This  is 
consistent  with  a  change  in  lidar  extinction  to  backscatter  ratio  noted  by  Platt  and 
Dilley  (1981).  We  have  exploited  this  dependence  to  simplify  our  treatment  of 
the  temperature  dependence  of  the  cirrus  backscatter,  and,  hence,  the  assignment 
of  simulated  wind  errors. 
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To  calculate  backscaltcr  coefficients  we  adopt  the  optical  constants  for  ice 
from  Warren  (1984).  We  have  evaluated  the  backscattcr  and  extinction  coef¬ 
ficients  corresponding  to  the  temperature  dependent  size  distributions  described 
above  using  the  Mie  theory  algorithm  of  Shcttle  (pers.  comm.)  Small  panicle 
contributions  were  added  to  complete  the  size  distribution  data.  Although  Mie 
theory  is  strictly  applicable  to  spherical  particles,  a  comparison  of  asymmetry  fac¬ 
tors  calculated  by  Takano  and  Lion  (1989)  for  hexagonal  ciA'stals  and  equivalent 
ice  spheres  suggests  that,  for  our  purposes,  the  error  is  not  significant.  Figures 
7  and  8  illustrate  the  extinction  (km  ’)  and  backscattcr  (km’’  sr‘*),  respectively, 
which  were  used  in  the  lidar  sounder  simulation. 

Modeled  cirrus  back.scatter  values  for  "cold"  and  "wann"  cirrus  at  the  lidar 
wavelength  of  2.1  /zm  correspond  to  about  1.5  x  10'^  and  5.0  x  10'^m'*sr'\ 
respectively.  The  appropriate  visible  values  are  approximately  three  times  larger, 
corresponding  to  "weak"  and  "medium"  cirrus.  Weak  cirrus  may  be  visible  with 
sufficient  vertical  thickness,  while  medium  cirrus  will  usually  be  visible. 

These  cirrus  back.scatter  parameters  as  well  as  the  appropriate  marine  aerosol 
and  liquid  water  cloud  models  were  provided  as  input  to  the  lidar  simulation 
model  at  Gli  AstroSpace  Division  (D.  Hogan,  pers.  comm.)  who  performed  the 
line  of  sight  and  horizontal  wind  vcl(x:ity  error  calculations.  Lidar  wind  errors 
were  then  evaluated  based  on  the  nominal  system  characteristics  listed  in  Table  6. 


Tabli'  (>.  I,i<I.ir  .SvMcin 


Faramctcr 

Value 

Wavelength 

2.1  fim 

Sensor  altitude 

824  km 

Scan  angle 

45  degrees 

Scan  pericxl 

S  s 

Pulse  energy 

1  J 

Pulse  duration 

.7  ,zs 

Noise  bandwidth 

2.2  MHz 

Pulse  rate 

10  Hz 

Optic-, ,  diameter 

.5  m 

Total  efficiency 

0.75 
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To  obtain  error  estimates  loi  the  total  measured  horizontal  wind,  shot  aver¬ 
aging  within  a  100  km  square  area  with  1  km  vertical  resolution  was  performed, 
and  typical  profiles  of  temperature  and  relative  humidity  were  assumed.  Sce¬ 
narios  considered  include  clear  skies  (boundary  layer  winds  only),  water  cloud 
only  (cloud  top  winds  only),  cirrus  clouds  (upper  level  winds  and  boundary  layer 
winds  with  reduced  accuracy),  and  cirrus  cloud  over  lower  clouds  (upper  level 
winds  and  cloud  top  winds  with  reduced  accuracy).  We  found  that  for  the  range 
of  backscatter  considered,  in  no  ca.se  did  the  estimates  of  horizontal  wind  velocity 
error  exceed  1  m.s"',  and  often  tliey  were  less  than  this.  Therefore,  we  adopted  a 
conservative  approach  and  assigned  a  1  ms“'  error  to  all  cloud/acrosol  scenarios. 
However,  note  that  this  value  is  repre.sentative  of  the  best  portion  of  the  satellite 
swath  whose  total  width  is  approximately  1800  km  (i.e.  halfway  between  the 
subtrack  and  the  swath  edge);  the  actual  horizontal  wind  measurement  error  will 
increase  near  both  the  swatli  edge  and  the  subtrack,  which  we  discuss  below. 


5.3.2  Measurement  Error  as  a  Function  of  Distance  to  Subtrack 


Recall  that  the  simulated  lidar  wind  measurements  assume  a  conical  scan,  as 
in  Figure  9  ,  with  a  radius  r.  The  lidar  only  measures  wind  in  the  radial  direction 
and  two  shot  averaged  measurements  are  needed,  one  in  the  forward  direction  Vf 
and  one  in  the  backward  direction  Vh,  to  define  the  total  wind  vector  V  Vj^). 
We  assume  throughout  that  all  measurement  errors  are  unbiased  and  random  wiUi 
no  correlations  in  the  vertical  or  horizontal.  If  the  measurement  error  on  the  radial 
velocity  is  6Vr,  we  wish  to  find  the  resulting  error  in  V  (5V0,  and  its  dependence 
on  the  distance  of  the  measurement  to  the  subtrack  line  (d). 

Projection  of  the  true  wind  onto  the  forward  and  backward  mea¬ 

surement  direction  results  in  the  following  transformation: 


I'V  -■  Fji  cosQ  -i  Vj  .sin  a 

.  .  .  (.5.3) 

I'i  -  V’ji  cos  a  -j-  sin  a 

Solving  for  Iji  and  V,  and  assuming  6\'f  --  SV'i,  -  6Vr,  the  errors  on  the 
components  of  V'  arc  then: 


6l’,i  .-rr  6Vv/^\/2 cos 
6  F I  --  6Vr/  y/2  si  n 


The  total  error  is  then 


(5.4) 


2S 


SV  =  +  8Vl  =  8Vrl(y2  sin  a  cos 

r2  (5.5) 

=  -^====:--=/)l4 

-  d2) 

The  minimum  is  6V  =  \/28Vr  at  d  —  r/\/2,  i.e.  at  45  degrees.  Figure  10 
shows  the  error  amplification  factor  8Vl8Vr  in  terms  of  dir. 

We  apply  this  analysis  to  the  specific  configuration  used  in  this  experiment 
to  obtain  values  of  the  relevant  parameters.  Recall  that  the  1  ms"’  error  value 
previously  derived  from  the  lidar  simulation  model  was  for  the  total  wind  error 
at  the  best  point  in  the  satellite  swath.  This  corresponds  to  an  error  in  the  radial 
velocity  measurement  of  .707  ms"’ .  Table  7  contains  the  appropriate  parameters 
for  the  simple  lidar  experiment  as  a  function  of  distance  to  the  satellite  subtrack 
(or  alternatively,  n,  the  scan  position  index).  Amp  is  the  error  amplification  factor 
for  the  total  measured  wind  vector,  given  by  Amp  =  (\/2sin  a  cos  a)  Note 
that  nadir  .scan  data  ju'e  included  but  not  used  during  the  actual  data  simulation. 

Tabic  7.  Aiiiplificaiion  in  Radial  Wind  Error  vs.  Disiancc  to  Subtrack  for  current  OSSE. 


Scan 

Position(n) 

d(km) 

d/r 

ofdcgrees) 

Error 

Amplification 

1 

-628.3 

-0.70 

-44.3 

-1.41 

2 

-375.8 

-0.42 

-24.7 

-1.86 

3 

-177.9 

-0.20 

-11.4 

-3.65 

4 

0. 

0. 

0. 

- 

5 

177.9 

11.4 

3.65 

6 

375.8 

■M 

24.7 

1.86 

7 

628.3 

0.70 

44.3 

1.41 

As  the  table  indicates,  the  SSM/T  scan  positions  correspond  to  data  only  from 
the  inner  half  of  the  satclliie  swath,  with  the  largest  measurement  errors  occurring 
for  measurements  directly  on  either  side  of  the  subtrack  line. 

Tlie  amplification  factor  shown  above  affects  the  measured  components  of  the 
wind  differently  at  any  given  point  in  time  depending  on  the  current  orientation 
of  the  satellite  subirack  with  resirect  to  lines  of  latitude  and  longitude,  'flie 
magnitude  of  the  errors  in  the  measured  along  and  cross  track  wind  components 
(ctjii.  5.4)  depeiui  only  on  scan  |H)sition.  However,  the  along  and  cross  track 
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direction  vectors  arc  constantly  updated  througliout  the  satellite  orbit  to  account 
for  the  changing  satellite  path.  For  example,  at  lower  latitudes  the  s  '.tellite  cross 
track  direction  will  be  nearly  parallel  to  latitude  lines,  while  near  the  poles  the 
cross  track  direction  will  be  roughly  parallel  to  meridians.  Thus,  random  errors  of 
the  proper  size  for  Vji  and  Vx  are  generated  and  then  are  transfonned  to  u  and  v 
measurement  errors  which  are  added  to  the  exact  components  sampled  previously 
from  the  nature  run. 

6.  IMPACT  RESULTS 

Evaluation  of  the  OSSH  impact  results  is  ba.sed  on  examination  of  synoptic 
maps,  the  time  evolution  of  analysis  and  forecast  error,  the  vertical  distribution  of 
analysis  and  forecast  error,  and  averaging  of  forecast  error  for  forecasts  initialized 
at  days  3,  5,  and  7  of  the  assimilation  period.  The  meteorological  fields  which  are 
u.sed  in  the  evaluation  include  the  500  liFa  geopotcniial  height,  850  hPa  vector 
and  zonal  winds,  and  the  850  hPa  relative  humidity.  Finally,  the  scaled  down 
lidtu"  results  are  compared  to  our  previous  W1ND.SA1'  results. 

6.1  Synoptic  Maps 

6.1.1  500  hPa  Geopotcniial  Height  Analysis  and  Forecast 

Figures  12  and  13  contain  the  height  field  analysis  and  analysis  error,  and 
the  48  h  forecast  and  forecast  error  at  500  liPa  for  both  the  CONTROL  and 
LIDAR  experiments  valid  at  0000  UTC  23  November  1979.  For  comparison,  the 
corresponding  Nature  (ields  are  shown  in  figure  11. 

In  the  northern  hemisphere  the  analy.ses  arc  nearly  identical,  with  little  impact 
on  the  dominant  wave  number  4  pattern.  Error  magnitudes  are  generally  less  than 
40  m  for  both  e>  ’'^  imenis.  One  location  where  the  lidar  data  seems  to  have  had 
an  impact  is  ne“  ’  ‘  0°W  where  the  negative  errors  of  "80  m  in  CONTROL  are 
reduced  to  less  .  ■  n  40  in. 

The  southein  hemisphere  sliows  a  much  stronger  impact  of  the  lidar  data,  with 
most  impact  seen  south  of  50°  .S  where  the  jet  stream  is  located.  For  example, 
the  broad  trough  between  9{FW  and  150°\V  is  better  analyzed  in  LIDAR.  In  this 
region  positive  errors  are  reduced  from  approximately  120-160  m  to  40-80  m, 
while  negative  errors  are  reduced  from  about  120  m  to  40  m.  Similarly,  the 
closed  low  near  90°E,  60°S  which  is  analyzed  160  m  t(x>  high  and  that  near 
30°E,  60°S  which  is  too  high  by  120  ni  in  CONTROL  are  both  better  analyzed 
in  LIDAR,  with  errors  reduced  to  40  m  and  80  m,  respectively. 
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Figure  12  5lK)  III’.!  ge<i|>ilenlial  lioiglil  .iti.ily'.i’:  valid  (K)  IJTC  23  November  1979.  (a).(b) 
CONTKf)L  analysis  anil  analysis  error.  (e),(d)  I.IDAR  analysis  and  analysis  error  for  (he 
norlhern  hemisplieio.  Conosiionding  soiitliem  liennspherc  plots  are  in  (e)-(h).  respectively. 
Contour  interval  is  80  in  lor  analyses,  and  >10  in  for  analysis  cnors.  (Continued  .  .  .  ) 
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Figure  12.  500  hi’a  gcopiuciiluil  lioighi  analysis  valid  (K)  UTC  23  November  1979.  (a).(b) 
CONTROL  analj.sis  and  analysis  error,  (c),(d)  LIOAK  analysis  and  analysis  error  for  the 
northern  hemisphere.  Corresponding  southern  hcmisplicre  plots  arc  in  (c)-(h).  respectively. 
Contour  interval  is  80  ni  for  analyses,  and  40  m  for  analysis  errors.  (Continued  . . .  ) 
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Figure  12.  500  hPa  gcopotcniial  height  analy.sis  valid  00  UTC  23  Novembtr  1979. 
(a),(b)  CONTROL  analy.<iis  and  analysis  error,  (c).(d)  LIDAR  analysis  and  analysis 
error  for  the  northern  hemisphere.  Corresponding  souUicm  hemisphere  plots  are  in 
(e)-(h).  respectively.  Contour  interval  is  80  m  for  analyses,  and  40  m  for  analysis  errors. 
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Figure  13.  5(K)  hl’a  48  hour  gt-o|>ok-nlij|  hcighl  aiulyMS  valiil  (M)  IPl'C  23  Novcmlicr  1979. 
(a),(b)  CONTRf)!,  analysiN  and  analysis  error,  (c).(d;  I.IDAK  analy.sis  and  analysis  error  for  llic 
northern  hemisphere.  Corresixindiiig  southern  hcmisplierc  plots  are  in  (c)  (li),  respectively. 
Contour  interval  is  80  m  for  analyses,  and  40  m  for  analysis  errors.  (Continued  .  .  .  ) 


Figure  13.  5()()  lil’a  48  liour  gciipolcnlial  height  analy.sis  valid  ()()  UTC  23  November  1979. 
(a),(b)  CONTROL  analysis  and  analysi.s  enor.  (c).(d)  LIDAR  analysis  and  analysis  error  for  the 
northern  hemisphere.  Corresponding  southern  hcniisplicrc  plots  arc  in  (c)-(h),  respectively. 
Contour  interval  is  80  m  for  analysc.s,  and  40  m  for  analysis  errors.  (Continued  .  . .  ) 
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Figure  13.  500  hl’a  48  hour  gcopotcntial  height  analysis  valid  00  UTC  23  November 
1979.  (a),(b)  CONTROL  analysis  and  analysis  error,  (c),(d)  LIDAR  analysis  and  analysis 
error  for  the  northern  hemisphere.  Corresponding  southern  hemisphere  plots  are  in 
(c).(h),  respectively.  Contour  interval  is  80  m  for  analyses,  and  40  m  for  analysis  errors. 


Tlie  forecast  fields  valid  at  the  same  time  show  tliat  the  DWL  data  has  litde 
impact  in  the  northern  hemisphere,  which  reflects  the  fact  that  both  forecasts 
started  with  very  similar  initial  states. 

However  in  the  southern  hemisphere  there  is  a  very  noticeable  impact  pole- 
ward  of  60°S.  In  general,  errors  are  reduced  in  magnitude  from  about  120-160  m 
in  CONTROL  to  roughly  40-80  m  in  LIDAR.  It  is  also  noteworthy  that  addition 
of  the  lidar  data  has  not  improved  phase  errors  in  the  forecast  since  the  same 
positive/negative  pattern  of  errors  is  retained  in  both  experiments.  For  example, 
the  trough/ridge  sequence  .south  and  east  of  Australia  in  Nature  is  moved  too 
slowly  in  both  CONTROL  and  LIDAR  forecasts  resulting  in  a  similarly  shaped 
error  pattern.  Since  the  same  model  is  used  in  both  forecast  runs,  errors  resulting 
from  inaccuracies  in  the  model  physical  parameterizations  would  also  be  present 
in  both  forecasts. 

6.1.2  850  hPa  Vector  Wind  Analy.sis  and  Foreca.sl 

Figure  14  contains  the  850  hPa  horizontal  wind  speed  analysis  and  vector 
wind  error  magnitudes  for  analysis  and  48  h  forecast,  while  figure  15  contains 
the  actual  horizontal  wind  vector  and  wind  vector  error  fields  for  analyses  and 
forecasts  valid  at  the  same  time,  00  UTC  23  November. 

For  the  analyses  most  of  the  data  impact  is  seen  in  the  southern  hemisphere, 
but  one  area  east  of  Newfoundland  associated  with  a  trough  over  eastern  North 
America  also  shows  an  improvement.  In  this  region  the  vector  wind  error 
maximum  magnitude  is  reduced  from  greater  than  30  m.s“'  to  approximately 
20  m.s"*. 

In  the  southern  hemisphere  the  impact  is  more  obvious.  Large  improvements 
are  seen  in  two  large  areas,  one  over  the  South  Pacific  Ocean  from  the  dateline  to 
South  America  near  60°S,  the  other  over  the  southern  Indian  Ocean  between  30°E 
and  90°E.  In  the  Pacific  much  of  this  improvement  is  due  to  the  closed  cyclonic 
circulation  near  60°S,  which  is  much  better  analyzed  in  LIDAR.  In  the 

Indian  Ocean  an  area  of  strong  westerly  winds  at  850  hPa  with  some  cyclonic 
curvature  is  poorly  analyzeil  in  CONTROL  but  much  improved  in  LIDAR.  In 
both  instances  areas  of  vector  wind  errors  of  20-30  ms“'  are  often  reduced  to  less 
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Figure  14.  850  hPa  global  liori/onial  wind  speed  analysis  and  vector  wind  error  magnitude  fields  valid  00 
UTC  23  November  1979.  (a),(b)  CONTROL  analysis  and  analysis  error,  (c).(d)  LIOAR  analysis  and 
analysis  error.  Corresponding  48  It  forecast  maps  valid  at  the  same  time  are  in  (e)-(h).  Contour  interval  is 
5  ms  '  starting  at  10  ms  '  lor  wiiul  speed,  and  starting  at  5  ms  '  for  the  error  m.ips.  (Continued  .  .  .  ) 
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I'ljiurc  14.  XSO  hl’a  global  liori/oni.il  wmil  .spcal  aiialy.sis  and  vector  wind  error  magnitude  fields  valid  00 
UTC  23  November  1979.  (a).(b)  CONTROL  analysis  and  analy.sis  error,  (c),(d)  LIDAR  analysis  and 
analysis  error.  Corresponding  48  h  forecast  maps  valid  at  the  .same  time  are  in  (c)-(li).  Contour  interval  is 
i  nis"'  starting  at  10  ms  ‘  for  wind  speed,  and  starting  at  .S  ms  ’  for  the  error  maps.  (Continued  .  .  .  ) 
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UTC  23  Novc! 
analysis  cnor.  O 
5  ms  *  starting  i 


1  igurc  14  hl’ii  j;1i)ImI  wiml  .spcotl  .iii.>I}Mn  jikI  \o<.-lor  wmil  oiror  ni.ip.niludc  ficIiK  valid 

{H)  UTC  2^  NovcinIxT  I'Pd  (.u.d))  CONI  KOI.  iiiulvMs  .ind  aii.dysiv  viior.  (c),(d)  I.IDAR  analysis 
and  analysis  otinr  C'virTispdndin);  4.S  li  loics.ist  maps  valid  ai  ihc  same  linn;  arc  'ii  (c)  (li)  Conlour 
interval  is  5  ms  '  si.ininj;  al  It)  ms  '  l.ir  svind  .speed,  and  Marling  al  >  ms  ‘  !oi  ilie  error  maps 
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f'lgurc  15  X50  lil\i  gloh.il  liori/oiii.ii  vmiuI  vector  .<iiat)sis  .ind  analysis  error  valid  (K)  ITFC  23 
November  1979.  (a),(b)  CONTROL  analssis  and  analysis  error,  (c).(d)  LIOAR  analysis  and  analysis 
error.  Concsponding  48  li  foreciisi  maps  valid  at  the  same  lime  arc  in  (c)-(h).  (Continued  .  .  ) 
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Figure  15.  850  hi’a  glob.il  hnri7oiUal  wind  vcelor  analy'!i<:  and  analysis  error  valid  00  UTC  23 
November  1979  (a),(b)  CONTROL  analy.sis  and  analysis  enor,  (c),(d)  LIDAR  analysis  and  analysis 
error.  Corresponding  48  li  forecast  maps  valid  at  llic  same  time  are  in  (e)-(li).  (Continued  .  .  .  ) 
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Fifjiiic  IS  KSd  hi’.i  j;li)lul  lioii/oiil.i!  wm,!  vn  loi  .lIl.l!yM^  .iml  .iiMly-is  crioi  v.ilui  (K»  UTC 
2?  NovimiiIh.t  !•>/<»  (.0.(I>1  I'ONIKOI.  .m.ilyviv  .iiul  .iii.ilysis  crioi.  (-,  ).(<1)  l.lDAk  .iii.iIvms 
and  jii.i!yM\  ciioi  ('Diti-si'DiKlitii;  IX  li  iii.ips  valid  al  die  same  liiiio  aie  ill  le)  (li). 


than  10  ins"'.  Forecast  maps  show  a  similar  pattern  with  impact  concentrated 
in  the  same  southern  hemisphere  locations.  Interestingly,  the  CONTROL  48  h 
forecast  error  seen  over  the  Indian  Ocean  is  somewhat  smaller  in  magnitude  than 
that  seen  in  the  analysis  valid  at  the  same  time.  This  is  consistent  with  the  fact  that 
compared  to  other  LIDAR  850  hPa  wind  forecasts,  the  error  growth  rate  for  the 
forecast  initialized  on  21  November  is  anomalously  small,  while  coincidentally, 
the  analysis  error  was  increasing  from  21  to  23  November  (sec  figure  16). 

6.1.3  850  hPa  Relative  Humidity  Analy.sis  and  Forecast  Error 

Figure  16  contains  the  relative  humidity  analysis  and  48  h  forecast  errors  for 
both  CONTROL  and  LIDAR  again  valid  for  (K)  UTC  23  November  1979.  As  in 
the  geopotential  height  fields,  nearly  all  impact  is  seen  in  the  southern  hemisphere 
extratropics,  south  of  3()°S.  In  particular,  large  areas  of  negative  error  in  the  south 
Pacific  and  Atlantic  with  levels  greater  than  50  percent  are  reduced  in  LIDAR 
analysis.  In  both  experiments  there  .seems  to  be  a  bias  toward  negative  errors 
rather  than  positive.  This  may  be  a  result  of  surface  moisture  parameterizations 
in  the  GSM  physics  package. 

In  the  forecast  maps,  impact  is  also  confined  to  the  southern  hemisphere  with 
large  improvement  seen  over  the  southern  oceans  and  Antarctica.  These  results  for 
the  humidity  analyses  and  forecasts  are  similar  to  those  obtained  in  the  WINDSAT 
OSSE  and  are  consistent  with  the  idea  that  improved  specification  (or  prediction) 
of  the  horizontal  wind  field  leads  to  a  better  analysis  (or  forecast)  of  a  passive 
tracer  like  water  vapor.  In  cs.sencc,  the  model  humidity  field  quickly  adjusts  to 
the  fields  of  mass  and  motion. 

6.2  Analy.si.s  and  Forcca.it  Krror  as  a  Function  of  Time 

6.2.1  500  hPa  (icopotcnlial  Height  rnis  Error 

Figure  17  shows  5(X)  hPa  mis  height  errors  of  the  7  day  analyses  and 
forecasts  for  several  experiments.  Errors  are  shown  for  the  globe  as  well  as 
the  northern  hemisphere,  southern  hemisphere  and  trojiical  subregions.  The 
experiments  plotted  are  CON'fROL  and  LIDAR,  and  for  comparative  purposes, 
the  previously  conducted  STALSAf  (control)  and  WINDSAT  OSSEs. 


Higyro  16.  K.SO  lil’.i  global  rcl.ilivo  liuir.uiiiy  analyMN  ciior  aiul  -IK  li  loroca\l  cnoi  valid  (X)  UTC  23 
No\cml)cr  1979.  (^).(h)  C‘().N’i'IU7l.  aiMlyM.-.  and  foro(..i.N!  error,  H).fd)  LIDAR  ,-uialy.Ms  error  and  forecast 
error  Carnlour  inle:\  al  is  2^  ]is'ks'ii1  w  illi  /s-io  line  oiiiilU’d  fsegaiise  values  are  d.tslied.  (Coiilinued  .  .  .  ) 


Figure  16.  850  lil’a  glolul  relative  huiiiuiily  aiialy.'ii';  error  ami  -18  li  forecast  error  valid  00  UTC  23 
Novcmlier  1970.  (a),(b)  CONTROL  analy.sis  and  foreca.st  error,  (c),(d)  LIDAR  analysis  error  and 
forecast  error.  Contour  interval  is  25  percent  with  /eio  line  omitted  Negative  values  are  dashed 
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Time  (Julian  day) 


b. 


Time  (Julian  day) 


Figure  17.  5(H)  hi’a  gcopolcnlial  liciglil  analysis  and  forceasi  errors  averaged  lu)ri/.ontaIly.  (a) 
Cilobal.  (b)  Northern  heinisplicrc,  (c)  Tropics,  (d)  Southern  hemisphere.  Shown  are  CONTROL 
(dot),  LIDAK  (solid).  STATSAT  (dash),  and  WINDSAT  (dash-dot -dot -dot).  (Continued  . ..  ) 
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With  the  exception  ol‘  the  tropics,  where  heiglit  and  wind  fields  are  not 
geostrophically  coupled,  we  see  that,  relative  to  CONTROL,  LIDAR  reduces 
height  analysis  errors  substantially.  This  is  most  dramatic  in  the  .southern  hemi¬ 
sphere  where  the  height  error  is  reduced  from  approximately  45  m  to  20-25  m. 
In  the  northern  hemisphere  where  the  conventional  observation  network  is  more 
extensive  the  improvement  is  smaller  -  from  15-20  m  to  about  12  m.  Globally 
the  reduction  is  from  25-30  m  to  roughly  12  m. 

Forecast  errors  are  also  reduced.  In  the  southern  hemisphere  extratropics 
the  accuracy  of  the  day  3  forecast  has  been  improved  relative  to  CONTROL  by 
an  equivalent  increase  in  forecast  range  of  36  hours.  Globally,  an  increase  in 
forecast  accuracy  relative  to  CONTROL  of  1 2  hours  is  seen,  however  northern 
hemisphere  forecasts  are  only  minimally  affected. 

Comparing  analysis  errors  for  the  LIDAR  and  WINDSAT  experiments  we 
see  that  depending  on  the  region,  WINDSAT  errors  are  either  smaller  or  larger 
than  LIDAR.  Since  each  OSSL  was  run  with  a  different  control  experiment  we 
can  compare  the  improvement  of  each  experiment  relative  to  their  respective 
control  runs.  This  indicates  that  WINDSAT  generally  shows  a  larger  relative 
improvement  over  STATSAT  than  docs  LIDAR  vs.  CONTROL.  For  example, 
the  southern  hemisphere  results  show  a  reduction  from  approximately  45  m  to 
12-15  m  with  WINDSAT,  compared  with  the  45  m  to  25  m  reduction  noted 
above.  In  the  northern  hemisphere  the  WINDSAT  improvement  over  STATSAT 
is  about  15  m  compared  to  only  5  m  in  LIDAR.  We  would  expect  a  larger  impact 
with  WINDSAT  data  since  the  instrument  was  assumed  to  return  more  complete 
measurements  throughout  the  clear  troposphere.  In  addition,  the  CONTROL 
analysis  is  generally  improved  reiative  to  STATSAT,  so  there  is  less  room  for 
improvement  due  to  the  additional  lidar  wind  data. 

6.2.2  850  hPa  Zonal  Wind  rms  Ern»( 

Figure  1 8  contains  the  analysis  error  of  tlie  850  hPa  u  component  wind  as  a 
function  of  time  during  the  assimilation  period.  Globally  the  error  is  reduced  from 
approximately  4-4.5  ms"'  to  3-3.5  ms"’  rms.  In  terms  of  region,  the  northern 
hemisphere  shows  only  veiy  small  improvement  over  CONTROL,  less  than  .5 
ms"’  to  -3  ms"’  rms.  In  the  tropics  the  reduction  in  error  is  also  small  but 
increases  through  the  assimilation  period  so  that  by  days  6  and  7  (Julian  days 
328  and  329),  the  rms  error  goes  from  about  4  ms"’  to  3.5  ms"’.  Apparently 
the  tropical  wind  variance  is  so  small  that  there  is  little  room  for  improvement. 
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Figure  18.  850  lil’a  ii  cdiniioncm  wind  analy.<iis  and  foreca.<;t  crroi.s  averaged  liori/ontally.  (a) 
Global,  (b)  Noriliern  lieinisphcrc,  (c)  Tropics,  (d)  Southern  hcniisplierc.  Shown  arc  CONTROL 
(dot),  LII3AR  (.solid).  STATSAT  (dash),  and  WINDSAT  (da.sh-doidol-dot).  (Continued  ...  ) 
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l-'iguro  18.  850  lil’a  ii  loinpoiu'iil  wind  iln.lly^l^  .iiul  l<in.-t.isl  cirois  .iveragcd  hori/onlally. 
(a)  Cilobal,  (b)  Nortliorii  ln;inispln;ic,  (c)  Tropics,  (d)  Soulliuni  hcmispbcro.  Shown  arc 
CONTROL  (dol),  LIDAU  (.solid).  STATSAT  (d.i.sli).  and  WINIXSAT  (dash-dol-dot-do!). 


even  with  relatively  high  quality  wind  data.  In  the  southern  hemisphere  we  again 
find  the  lidar  data  impact  to  he  much  greater.  In  CONTROL  nns  errors  arc 
4.5  and  5.5  ms  ’  at  the  start  and  end  of  the  assimilation  ixjritxl,  respectively. 
But  in  LIDAR  these  are  reduced  to  approximately  .'f.5  and  3  ms"’,  respectively. 
Forecast  accuracy  is  also  seen  to  improve,  particularly  in  the  southern  hemisphere; 
increases  in  predictability  arc  as  large  as  48  h. 

Compared  to  WlNDSAl',  improvements  are  mt't :  modest,  generally  no  greater 
than  one  half  the  magnitude  of  the  improvements  seen  in  WINDSAT  vs.  STAT- 
SAT.  Again,  the  comparison  is  somewhat  ambiguous  since  (1)  the  CDW  data  were 
modified  to  induce  horizontal  correlations,  and  (2)  the  01  had  been  significantly 
altered. 

6.2.3  2(M)  hPa  Zonal  Wind  rm.s  Krror 

Errors  for  the  200  hPa  u  component  wind  are  in  figure  19.  Impacts  here 
generally  resemble  those  at  850  hPa.  That  is,  almost  no  impact  in  the  northern 
hemisphere,  small  impact  in  the  tropics  (<  .5  m/s  improvement),  and  larger 
impact  in  the  southern  hemisphere.  Here  the  error  is  reduced  from  about  6  m/s 
rms  to  approximately  4-5  m/s.  Increases  in  predictability  with  the  DWL  data 
are  somewhat  smaller  than  at  850  hPa  with  improvements  in  the  range  12  to 
24  h.  Finally,  the  LIDAR  relative  impacts  are  generally  smaller  than  those  for 
WINDSAT. 

6.2.4  850  hPa  Relative  Humidity  rms  Error 

Globally,  improvement  over  CONTROL  in  the  analysis  error  for  relative 
humidity  (figure  20)  is  small.  Since  both  the  northern  hemisphere  and  tropics  show 
essentially  no  impact  most  of  this  improvement  is  concentrated  in  the  southern 
hemisphere.  Here  there  is  a  positive  impact  which  varies  with  time  from  less  than 
1  percent  at  the  start  of  the  assimilation  pcritxl,  to  roughly  3  percent  by  day  7. 
Forecasts  in  the  southern  hemisplu  re  are  improved  with  predictability  increases 
of  24  to  72  hours. 

6.3  Analysis  and  Forecast  Error  as  a  Function  of  Pressure 

Analysis  errors  were  also  averaged  over  latitude,  longitude,  and  over  the 
last  5  days  of  the  assimilation  period  to  obtain  vertical  profiles  of  error  for  the 
northern  hemisphere,  .southern  hemisphere,  and  tropical  subregions.  Plots  for 
geopotential  height,  vector  wind,  relative  humidity,  and  temperature  are  found  in 
figures  21-24,  respectively. 
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Figure  19.  200  hl’a  u  coniponoiii  wind  ainily.sis  and  forecisl  errors  averaged  horizontally,  (a) 
Global,  (b)  Northern  hemisphere,  (c)  Tropies,  (d)  .Southern  hemisphere  Slmwn  are  CUNTROL 
(Jot).  LIDAR  (solid),  STAT.SAT  (dash),  .and  WIND.SAT  (da.sh  dot-dot  dot).  (Continued  .  . .  ) 
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Figure  19.  200  hi’a  u  coniponcnl  wind  analysis  and  forccasl  errors  averaged  horizontally, 
(a)  Global,  (b)  Northern  hemisphere,  (c)  Tropics,  (d)  Southern  hemisphere.  Shown  are 
CONTROL  (dot).  LIDAR  (.solid).  STATSAT  (dash),  and  WINOSAT  (dash-dot-dot.dot). 
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Figure  20.  850  hi’j  relative  luinudily  aiialysi.s  and  forecast  cnors  averdgc<l  horizontally,  (a) 
Global,  (b)  Northern  hcniisphere.  (c)  Tropics,  (d)  .Southern  heinispheie.  Shown  .ve  CONTROL 
(dot).  LIOAR  (solid).  STAT'SAT'  (dash),  .md  WINOSAT  (dash-dol-dot  dot).  (Continued  .  . .  ) 
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llic  assimilalinn  period,  (a)  OlobjI.  (h)  Notlhem  honii.splierc,  (i  )  Tropics,  (d)  Southern  hcmtsphcrc.  Shown  arc 


CONTROL  (dot).  LIDAK  (solid).  STATSAT  (da.sh).  and  WINIXSAT  (<lash  dot-<k>t-dol).  (Continued  . . .  ) 


PiKiirc  21.  Kcl.iiivc  luimidily  nnatysix  error  vs.  pressure  nveriiKcd  liori/oninlly  and  over  the  last  .5 
days  of  the  assimilation  period,  (a)  Northern  hemisphere,  (b)  Tropics,  (c)  Southern  hemisphere. 
Shown  lire  CONTROL  (dot).  LIDAR  (solid).  STATSAT  (dash),  and  WINDSAT  (dash-dot-dot-dot). 
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f-igurc  24.  Temperature  aii.ilysis  error  vs.  pressure  averaged  liori/.ontally  and  over  the  last  5  days  of 
the  assimilation  period,  (a)  riIob.il.  (b)  Northern  hemisphere,  (c)  Tropics,  (d)  Southern  hemisphere. 
Sliown  arc  CONTROL  (dot).  I.IDAR  (solid).  STATSAT  (d.ish).  and  WINDSAT  (dash-dol-dol-dol). 
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6.3.1  Gcopotcntial  Height 

As  seen  in  the  spatially  averaged  analysis  errors  of  section  6.2  the  largest 
impact  is  found  in  the  southern  hemisphere.  Here  errors  are  reduced  from 
about  40-50  m  in  CONTROL  to  roughly  30  m  between  1000  and  100  hPa. 
Compared  to  WINDSAT  the  impact  in  LIDAR  is  smaller,  even  accounting  for 
the  improved  analysis  system.  The  one  exception  to  this  is  above  100  hPa  where 
both  CONTROL  and  LIDAR  generally  have  smaller  errors  than  either  STATSAT 
or  WINDSAT. 

Vertical  profiles  of  the  height  analysis  errors  show  the  improvement  of  LIDAR 
over  CONTROL  to  be  relatively  constant  with  pressure.  Tliis  indicates  that  most 
of  the  improvement  to  the  height  analyses  is  barotropic  (i.e.  layer  thicknesses  are 
not  improved)  and  that  most  of  the  impact  is  a  reflection  of  an  improvement  in 
the  surface  pressure  analysis,  which  affects  all  levels  equally. 


6.3.2  Vector  Wind 

Similar  to  the  earlier  wind  analysis  results,  the  strongest  effect  of  the  DWL 
data  is  seen  in  the  southern  hemisphere  where  the  rms  vector  wind  error  is  reduced 
from  6.5-9.5  ms“*  to  about  5-9  ms“’,  depending  on  level.  In  other  regions  th  '. 
impact  is  no  greater  than  roughly  1  ms"*.  Again  it  appears  that  WINDSAT  impacts 
are  larger  with  mis  errors  falling  to  3-4.5  ms"*. 


6.3.3  Relative  Humidity 

These  plots  indicate  that  in  all  regions  of  the  globe  which  were  examined, 
specification  of  the  relative  humidity  field  is  improved  relative  to  CONTROL.  The 
smallest  positive  impact  is  .seen  in  the  northern  hemisphere  and  tropics  where  the 
LIDAR  analysis  error  is  generally  smaller  than  CONTROL  by  less  than  1  percent. 
In  t’.^  southern  hemisphere  the  observed  impact  is  much  larger.  At  and  above 
850  i.r'a  the  rms  analysis  error  is  reduced  2.5  to  3  percent,  while  at  1000  hPa  the 
reduction  falls  to  .5  percent.  Between  700  and  300  hPa  WINDSAT  still  appears 
to  have  a  much  larger  impact  in  all  regions  reducing  the  rms  error  to  about  17 
percent  in  the  northern  hemisphere  and  20-22  percent  in  the  .southern  hemisphere. 
These  results  are  consistent  with  those  seen  in  the  horizontally  averaged  statistics 
and  in  the  synoptic  plots  above. 
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6.3.4  Temperature 

Impacts  upon  the  trojx)spheric  temperature  analyses  are  small  for  all  experi¬ 
ments.  Below  150  liPa  most  differences  are  less  than  .5  K.  with  the  largest  impacts 
seen  in  the  southern  hemisphere.  However  as  seen  in  the  height  analysis  errors, 
both  CONTROL  and  LIDAR  are  superior  to  STATS  AT  and  WINDS  AT  above  150 
hPa,  with  typical  reductions  of  error  on  the  order  of  1  to  4  K.  This  might  be  a 
reflection  of  improvements  made  in  the  handling  of  raob  height  data  in  the  01, 
which  tend  to  affect  both  the  height  analyses  and  the  temperature  analyses. 

6.4  Averaged  Forecast  Errors 

Errors  associated  with  the  day  3,  5,  and  7  forecasts  were  combined  to  produce 
a  single  averaged  nns  forecast  error  curve  for  each  experiment  and  forecast 
variable.  Results  for  5(X)  hPa,  200  hPa  zonal  wind,  850  hPa  zonal  wind,  and 
850  hPa  relative  humidity  are  found  in  figures  26  to  29,  respectively. 

In  general,  the  averaged  errors  corroborate  the  features  of  individual  forecast 
errors  discussed  in  the  earlier  section;  forecast  impacts  are  large.st  in  the  southern 
hemisphere  and  negligible  in  the  northern  hemisphere.  In  the  case  of  the  500 
hPa  height  forecast  errors  the  averaged  improvement  of  LIDAR  vs.  CONTROL 
is  of  the  order  2-2.5  days  in  the  .southern  hemisphere  and  about  1  day  globally. 
WINDSAT  errors  are  better  still  by  1  day  in  the  southern  hemisphere  and  about 
12  h  globally.  Zonal  wind  foreca.st  errors  also  show  the  greatest  impact  of 
the  DWL  data  to  be  in  the  south.  Here  the  averaged  improvement  in  forecast 
accuracy  is  1-1.5  days  at  200  hPa  and  roughly  2  days  at  850  hPa.  Globally 
the  corresponding  improvements  are  .5  and  1  day,  re.spectively.  WINDSAT  still 
yields  superior  forec.ists,  improving  over  LIDAR  in  the  southern  hemisphere  by 
about  1.5  days.  Finally,  relative  humidity  forecasts  in  LIDAR  are  better  than 
CONTROL  by  1-1.5  days  in  the  .southern  hemisphere,  and  WINDSAT  errors, 
which  are  initially  lan>er  than  LIDAR  in  the  south  eventually  become  smaller 
than  LIDAR  by  approximately  12  h.  In  the  tropics  and  northern  hemisphere 
LIDAR  does  not  improve  forecast  errors  relative  to  CONTROL  and  WINDSAT 
errors  in  this  region  are  actually  larger  than  cither  CONTROL  or  LIDAR. 
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Figure  25.  Number  of  lidar  return.':  as 
UTC  19  November  1979.  A  sigma  vali 
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l  i};iirc  26.  .“ilM)  lil’.t  Ik'ij;IiI  (orcc.isl  i-tmr.s  ;ls  a  liincluin  <il  liim-  avi-iajjcil  over  llio  day 
3.  .i.  and  7  forecasts  (a)  Globat,  (b)  Norllierii  Hciiii.spliere,  (e)  Tropics,  (d)  Soutliorn 
ilcmisplicrc.  lixperiiiiems  shown  arc  l.IIMR  (solid).  CONTROL  (doi),  STATSAT  (da.sh),  and 
WINDSAT  (dash-dol-dol-dol).  Fic  27  through  29  have  the  same  fornial.  (Continitcd  . .  ) 
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Figure  26.  .“iiK)  lil’j  height  foreca.<;i  errors  as  a  function  of  time  averaged  over  the  day 
3,  and  7  forecasts,  (a)  Global,  (b)  Nortlicrn  Hemisphere,  (c)  Tropics,  (d)  Southern 
Hemisphere.  Experiments  shown  are  LIDAK  (solid).  CONTROL  (dot).  STATSAT 
(dash),  and  WINDSAT  (dash-dot-dot-dot).  Figs.  27  through  29  have  the  same  formal. 


Zonal  wind  (m/sec)  Zonal  wind  (m/sec) 


0.0  0.5  1.0  1.5  2.0  2.5  3.0  3.5  4.0 
Forecast  length  (days) 


Forecast  length  (days) 


p!"urc  27.  850  hi’a  /x>iial  wmk!  foiccast  errors  as  run<-lion  of  lime  avera.'cd  over 
day  3.  5,  and  7  forecasts.  Sec  fig,  26  for  formal  and  labelling.  (Conlinecd  .  . . 
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Figure  27.  850  lil’a  zonal 
over  llie  day  1,  5.  and  ' 
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Figure  28.  20()  liP.i  /on.il  wind  lorcca.-Jl  errors  its  a  funclion  of  lime  averaged  over  the 
day  3,  5,  and  7  forccasis  .See  lig  26  for  format  and  laltelling  (Continued  .  .  .  ) 
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f'lgurc  28.  2(H)  lil’a  zonal  wind  forcc.i.sl  errors  as  a  lunctinn  of  lime  averaged 
over  the  day  3,  5,  and  7  forcca.sls.  Sec  fig.  26  for  format  and  labelling. 
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Figure  2‘t.  K.'iO  111’.!  iclalivo  hiiiiiiclily  loroea.sl  errors  as  a  riitiction  ol  time  averaged  over 
the  (lay  0.  5.  and  7  forecasis.  See  fig.  ic.  lor  foimal  and  lalxilling.  (Continued  .  .  .  ) 
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6.5  Comparison  with  WINDSAT 

Because  WINDSAT  and  LIDAR  were  run  using  different  analysis  systems, 
gauging  relative  data  impact  is  rather  ditficult.  The  control  experiment  corre¬ 
sponding  to  WINDSAT  (i.e.  STATSAT)  had  substantially  larger  errors  in  the 
northern  hemisphere  than  the  control  experiment  using  the  new  analysis  system 
(CONTROL)  -  the  difference  in  northern  hemisphere  500  hPa  height  errors  be¬ 
tween  CONTROL  and  STATSAT  were  larger  than  the  WINDSAT  impact  using 
the  old  01.  Thus,  although  the  northern  hemisphere  LIDAR  impact  was  negligi¬ 
ble  and  the  WINDSAT  impact  was  large  and  positive,  no  definitive  conclusions 
can  be  drawn  about  the  relative  impact  of  the  two  lidar  .sounders  in  the  northern 
hemisphere,  since  the  WINDSAT  impact  is  certain  to  be  much  reduced  for  the 
new  analysis  system.  The  overestimate  of  the  WINDSAT  impact  in  the  northern 
hemisphere  is  due  to  the  fact  that  in  the  old  analysis  system  the  WINDSAT  data 
provided  the  only  infonnation  about  the  surface  pressure  gradient. 

In  the  tropics  and  southern  hemisphere  tlte  differences  between  CONTROL 
and  STATSAT  are  smaller,  allowing  a  more  meaningful  comparison  between  the 
two  impact  experiments.  The  zonal  wind  itnpact  in  the  tropics  was  essentially  zero 
for  LIDAR,  but  it  was  large  and  positive  in  WINDSAT.  The  new  analy.sis  system 
also  reduced  errors  compared  to  STATSAT,  but  by  much  less  than  WINDSAT, 
so  that  the  WINDSAT  impact  can  reasonably  be  expected  to  remain  positive 
even  for  the  improved  analysis.  In  the  southern  hemisphere,  both  the  height 
and  zonal  wind  errors  were  little  affected  by  the  new  analysis  system,  so  that 
WINDSAT  and  LIDAR  impacts  can  be  directly  compared.  In  terms  of  500  hPa 
heights,  the  reduction  in  analysis  error  in  LIDAR  is  l)etwcen  1/2  and  3/4  that  of 
WINDSAT,  and  the  LIDAR  increa.se  in  predictability  of  1.5  to  2  days  is  about 
1  day  less  than  that  of  WINDSAT.  The  zonal  wind  mises  at  850  hPa  show  a 
LIDAR  impact  that  is  half  the  magnitude  as  that  .seen  in  WINDSAT,  whereas  it 
is  slightly  smaller,  both  in  absolute  and  relative  terms,  at  2(X)  mb.  This  may  be  a 
direct  reflection  of  the  data  distribution  in  the  two  systems:  the  LIDAR  sounder 
only  has  data  at  200  hPa  if  there  is  cirrus  present  at  that  level,  whereas  marine 
boundary  layer  winds  were  assumed  available  everywhere  where  the  return  signal 
was  not  blocked  by  intervening  cloud  layers.  1-or  example,  figure  25  shows  the 
vertical  distribution  of  lidar  data  for  1200  I'TC  on  19  Novemlier  indicating  a 
sharp  peak  near  .sigma-.9,  and  much  fewer  observations  at  sigma=.2  -.3.  Since 
tliis  pattern  suggests  that,  to  some  extent,  differences  in  impact  between  LIDAR 
and  WINDSAT  are  a  function  of  data  availability  or  volume,  it  may  be  reasonable 
to  consider  the  use  of  two  simultaneous  DWL-bearing  satellites,  each  of  which 
can  fill  in  the  gaps  in  horizontal  coverage  left  by  the  other.  Relative  humidity 
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impacts  were  primarily  limited  to  the  southern  hemisphere  for  LIDAR,  unlike 
WINDSAT  which  also  improved  the  specification  in  the  northern  hemisphere. 
But  again,  WINDSAT  impacts  may  be  overstated  in  the  north  because  of  the 
analysis  sy.stcm  used.  It  is  likely  that,  run  with  the  updated  01,  WINDSAT  would 
yield  impacts  similar  to  LIDAR. 

7.  CONCLUSIONS 

We  have  conducted  an  OSSL  to  determine  the  potential  impact  of  a  simple  low 
energy  DWL  sensor  on  a  global  assimilation  and  forecast  system.  Measurements 
were  simulated  assuming  that  returns  would  mainly  be  available  at  or  near  cloud 
tops  and  in  the  marine  boundar>'  layer.  This  can  be  contrasted  with  measure¬ 
ments  from  LAWS  which  are  expected  to  be  more  complete  with  measurements 
throughout  the  clear  troposphere. 

We  have  found  that  measurements  from  such  a  scaled  down  DWL  have  a 
positive  impact  on  analyses  and  forecasts.  Analyses  errors  for  geopotential  height 
are  reduced  in  both  northern  and  southern  hemisphere  extratropics,  but  especially 
in  the  south  where  conventional  observations  are  few.  Analyses  of  horizontal  wind 
arc  improved  in  all  regions,  with  the  largest  impact  again  seen  in  the  southern 
hemisphere.  Improvements  in  forecasts  initialized  with  the  new  analyses  are  also 
found.  In  some  cases  rins  forecast  errors  of  height  are  improved  an  equivalent  of 
24-36  hours  in  forecast  time  relative  to  control  cases. 

The  major  results  of  the  OSSH  are  summarized  by  these  specific  findings: 

1.  Synoptic  maps: 

•  Reduction  in  analysis  and  forecast  errors  of  500  hPa  geopotential  height 
in  the  .southern  hemisphere  by  factors  of  2-3; 

•  Most  of  the  data  impact  on  vector  wind  analyses  and  forecasts  is  in  the 
southern  hemisphere  with  wind  vector  errors  of  20-30  m  s"'  reduced  to 
10  m  s"'; 

•  For  relative  humidity,  most  impact  also  in  the  southern  extratropics  with 
large  areas  of  errors  greater  than  50%  in  the  CONTROL  experiment 
reduced  in  the  LIDAR  analysis.  Similar  results  for  humidity  forecasts 
with  improved  specification  of  the  horizontal  motion  field  leading  to  a 
becer  moisture  forecast. 

2.  Time  evolution  of  forectist  and  analysis  error: 

•  Global  reduction  of  the  500  hPa  height  analysis  error  by  a  factor  of  two 
(from  25  m  to  12  m),  most  dramatic  in  the  southern  hemisphere.  A  global 
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increase  in  accuracy  for  a  day  3  forecast  by  the  equivalent  of  12  h  relative 
to  CONTROL  (36  h  in  the  southern  hemisphere  extratropics); 

•  A  thirty  percent  reduction  in  the  error  of  the  850  hPa  zonal  wind  com¬ 
ponent  with  a  forecast  predictability  improvement  in  the  southern  hemi¬ 
sphere  of  up  to  48  h; 

•  Consis»''nt  positive  impact  on  the  analysis  of  850  hPa  relative  humidity 
in  the  southern  hemisphere  with  forecast  predictability  increases  of  24  to 
72  h  compared  to  CONTROL. 

3.  Averaged  forecast  error: 

•  Global  improvement  of  the  500  hPa  height  forecast  (compared  to  CON¬ 
TROL',  by  1  day  (2-2.5  days  in  the  southern  hemisphere); 

•  Global  zonal  wind  forecast  error  improvements  of  0.5  and  1  day  at  850 
hPa  and  2(X)  hPa,  respectively.  (1-1.5  days  and  2  days,  respectively,  in 
the  southern  hemisphere); 

•  Improvements  in  the  850  hPa  humidity  forecasts  by  1-1.5  days. 

We  have  also  performed  an  OSSE  to  evaluate  the  impact  of  a  WINDSAT  type 
instrument  (Hoffman  ct  al.,  1990).  In  the  examination  of  southern  hemisphere 
results,  a  meaningful  comparison  of  WINDSAT  and  LIDAR  results  can  be  made. 

The  scaled  down  DWL  system  reduces  500  hPa  geopotential  analysis  errors  by 
1/2  to  3/4  that  of  the  WINDSAT  system.  Predictability  increases  for  LIDAR  are 
about  one  day  less  than  WINDSAT.  Zonal  wind  improvements  are  about  1/2  of 
WINDSAT.  However,  deficiencies  in  the  analysis  sy.stem  used  in  WINDSAT  tend 
to  exaggerate  the  positive  impact  of  a  WINDSAT  instrument  relative  to  LIDAR 
(sUiich  was  run  with  an  enhanced  analysis  system)  in  the  northern  hemisphere, 
making  intercomparison  in  tliis  region  more  difficult. 

To  a  large  degree,  results  obtained  here  substantiate  previous  data  impact 
studies  which  u.sed  simulated  DWL  data.  What  this  study  has  attempted  to  do  is 
show  that  a  low  energy  DWL,  which  may  be  viewed  as  a  lower  cost  alternative 
to  proposed  instruments  such  as  LAWS,  might  also  improve  the  global  analysis 
and  prediction  of  mass  and  motion  fields.  To  the  extent  that  the  results  we  have 
sliown  arc  valid,  such  a  scaled  down  DWL  should  receive  serious  consideration 
by  the  remote  .sensing  and  atmospheric  sciences  community  in  any  discus.sion  of 
future  spacc-based  observing  systems.  j 

'lliese  results  represent  a  lirst  attempt  to  assess  the  viability  of  a  simple  DWL 
concept  for  NWP.  As  siicii,  we  would  suggest  that  follow  up  studies  would  help 
c.xicnd  and  validate  the  results  sliown  here.  For  cxamjilc,  we  might  wi.sh  to 
determine  impact  sensitivity  to  various  orbital  inclinations.  Results  which  show  a 
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larger  impact  in  the  southern  hemisphere  would  seem  to  argue  for  a  more  highly 
inclined  orbit.  Additionally,  as  with  all  OSSEs,  the  issue  of  simulated  errors  is 
critical.  Our  measurements  assumed  unbia.scd  -I  m  s'’  nns  errors.  Further  off¬ 
line  studies  using  radiative  transfer  calculations  would  help  determine  whether  the 
level  of  accunicy  might  be  more  dependent  upon  cloud  microphysical  properties. 
Additional  consideration  is  also  needed  of  whether  some  cloud  and  aerosol  regimes 
might  lead  to  horizontally  correlated  errors. 
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